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Abstract 

The mean curvature flow describes the parabolic deformation of embedded branes in 
Riemannian geometry driven by their extrinsic mean curvature vector, which is typically 
associated to surface tension forces. It is the gradient flow of the area functional, and, 
as such, it is naturally identified with the boundary renormalization group equation of 
Dirichlet sigma models away from conformality, to lowest order in perturbation theory. 
D-branes appear as fixed points of this flow having conformally invariant boundary con- 
ditions. Simple running solutions include the paper-clip and the hair-pin (or grim-reaper) 
models on the plane, as well as scaling solutions associated to rational (p, q) closed curves 
and the decay of two intersecting lines. Stability analysis is performed in several cases 
while searching for transitions among different brane configurations. The combination 
of Ricci with the mean curvature flow is examined in detail together with several ex- 
plicit examples of deforming curves on curved backgrounds. Some general aspects of the 
mean curvature flow in higher dimensional ambient spaces are also discussed and obtain 
consistent truncations to lower dimensional systems. Selected physical applications are 
mentioned in the text, including tachyon condensation in open string theory and the 
resistive diffusion of force-free fields in magneto-hydrodynamics. 
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1 Introduction 



The theory of geometric flows is a modern subject of common interest in physics and 
mathematics. In abstract terms, one is led to consider continuous deformations of geo- 
metric quantities defined on Riemannian manifolds (of fixed topology) starting from some 
initial data and evolve it under a given set of parabolic equations with respect to the 
flow variable t, called time. When the driving term of the deformation is provided by the 
curvature, in various forms, the corresponding geometric evolutions are called curvature 
flows. These are naturally divided into two distinct classes, intrinsic and extrinsic curva- 
ture flows, depending on whether one deforms the metric on a Riemannian manifold M. 
by its Ricci curvature or a submanifold M embedded in M. by the associated extrinsic 
curvature. Since extrinsic geometry is more elementary than intrinsic geometry, as the 
simple example of planar curves illustrates, one expects that extrinsic curvature flows 
have longer history in science, as they do. The mathematical motivation for introducing 
geometric flows varies from one type to the other and the same also applies to their 
physical origin and diverse applications. However, what makes them worth studying is 
the undisputed fact that such a simple minded framework, based on geometric analysis, 
had far reaching consequences and led to ground breaking results in recent years. 

The main qualitative feature of curvature flows is their tendency to dissipate any pos- 
sible deviations from canonical geometries associated to fixed point configurations with 
special curvature. It is typical that intrinsic flows will deform the metrics towards con- 
stant curvature metrics, if they exist on a given Riemannian manifold. Likewise, extrinsic 
flows will deform the embedded submanifolds towards special configurations with pre- 
scribed extrinsic curvature, e.g., minimal submanifolds and generalizations thereof. This 
behavior is expected from parabolic equations on Riemannian manifolds, since they share 
common properties with the heat equation; actually, the latter is a local linear approxi- 
mation to the late stage evolution of geometries in the vicinity of the fixed points. These 
particular properties have turned geometric flows into a valuable tool for addressing a 
variety of long standing problems in differential geometry, such as the geometrization of 
manifolds in low dimensions, and many others. In practice, one should device suitable 
systems of flows for a given class of geometric data and follow their evolution towards 
configurations with prescribed curvature. In all cases there are many technical obstacles 
related to the possible formation of singularities along the flows, and their mathematical 
classification, which affect the long time behavior of solutions and need to be accounted 
for complete study. 
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The prime example of intrinsic curvature flow is the celebrated Ricci flow that de- 
forms metrics by their Ricci curvature tensor. It first arose in the physics literature as 
renormalization group equation for the coupling constant of non-linear sigma models de- 
fined in two dimensions, pQ. In this case, the deformation variable t is the logarithm of 
the world-sheet length scale of the field theory and the metric of the target space -viewed 
as generalized coupling- receives counter-terms that are computable by the perturbative 
renormalizability of the two-dimensional quantum field theory. Then, in this context, 
the Ricci flow describes the response of the target space metric to different energy (and 
hence length) scales of the quantum theory to lowest order in perturbation theory, [2J. 
When the renormalization group equation of sigma models is applied to spaces of positive 
curvature, it implies, in particular, that the quantum theory becomes asymptotically free 
in the ultra-violet regime, thus justifying the use of perturbative calculations at high 
energies. This is analogous to the asymptotic freedom exhibited by non-Abelian gauge 
theories in four space-time dimensions, [3], which, after all, motivated the study of renor- 
malization group flows in toy quantum field theories, such as non-linear sigma models in 
two dimensions. 

The derivation of the beta function equations for non-linear sigma models played a 
prominent role in the development of string theory, since fixed points of the Ricci flow are 
selected by the requirement of conformal invariance on the world-sheet. Generalizations 
in the presence of dilaton and anti-symmetric tensor fields were also considered in the 
literature, [U E], [6J, thus leading to coupled systems of beta function equations for the 
massless modes of closed strings which can be derived from an effective gravitational 
action in target space. Critical string theory in curved spaces is only concerned with 
the existence and construction of fixed point solutions to these equations. However, it 
was realized in recent years that the problem of tachyon condensation in closed string 
theory can be studied as off-shell process involving trajectories between different fixed 
points, via the renormalization group equations, in the weak gravitational regime. Thus, 
genuine running solutions of the Ricci flow, and its generalizations, are of great interest 
for exploring the problem of vacuum selection in closed string theory. 

The Ricci flow was introduced independently in the mathematics literature as new an- 
alytic tool to attack Poincare's conjecture and related geometric problems in three dimen- 
sions, [7j; in this context it also became known as Hamilton- Ricci flow. Since then there 
have been may important developments which are summarized in Ref. [H[9], together 
with the complete list of original contributions to the subject. Quite recently, Hamilton's 
programme was brought to completion by solving the long standing geometrization prob- 
lem of compact 3-manifolds by the Ricci flow, [ID] , in all generality; see also Ref. [TT] 
for an overview of this subject. The classification of singularities that may arise in the 
process and the introduction of appropriate entropy functionals for the Ricci flow have 
played important role in these studies. However, their relations to physics have not been 
entirely clarified so far. It also remains to understand in general mathematical terms the 
structure of the generalized Ricci flows, in the presence of anti-symmetric tensor fields, 
and their relevance to the formation of singularities. Finally, instanton corrections to the 
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beta function equations of non-linear sigma models, [12J, which are quite customary in 
physics and usually affect the nature of the infra-red fixed points, are still awaiting for 
their proper mathematical interpretation in the context of geometric analysis by Ricci 
flow. Non-trivial infra-red fixed points are known to exist by adding topological torsion 
terms with 9 = 7r in sigma models, [T3] . 

The prime example of extrinsic curvature flow is provided by the mean curvature flow 
of hypersurfaces that deform by their extrinsic mean curvature vector in the ambient 
space. It first arose in the physics literature as idealized (two-dimensional) model for the 
motion of grain boundaries in an annealing piece of metal, [HJ, by drawing analogies with 
the motion of soap bubbles and interfaces due to their surface tension; see also Ref. [15] 
for a similar model for the development of a surface groove by evaporation-condensation 
mechanism. The subject was subsequently generalized and put on firm mathematical 
base in Ref. [16]. Later, it flourished and became a whole area of intense study in math- 
ematics up to the present time; see, for instance, the recent textbooks [TFJ [181 tlH] and 
references therein. The mean curvature flow arises, in its original formulation, as gradient 
flow for the area functional of a hypersurface embedded in a fixed Riemannian manifold 
and as such it encompasses minimal surfaces among its critical points. Thus, it offers 
a new analytic tool in the framework of geometric analysis for studying minimal sub- 
manifolds, such as geodesies curves, and various isoperimetric problems associated with 
them. The structure of the singularities that may form on the way and the construction 
of entropy functionals for this flow are some of the main technical problems which are 
well investigated by now, as for the Ricci flow. 

Apart from its intrinsic mathematical interest, the mean curvature flow has several 
physical applications serving as local model for the evolution of interfaces and the den- 
dritic crystal growth, [20], the formation of labyrinthine patterns in ferro-fluids, [21J, the 
rendezvous problem for mobile autonomous robots, [22], and others that will be discussed 
in due course. There are also some variants of this flow, which will not be examined in 
this paper, that have led to ground breaking results in the mathematical theory of general 
relativity, [23]. However, what has been lacking, up to this day, was an account of the 
mean curvature flow in quantum field theory analogous to the Ricci flow. 

The primary aim of the present work is to describe in detail the field theoretic mani- 
festation of the mean curvature flow, and some of its generalizations, as boundary renor- 
malization group equations for Dirichlet sigma models defined on two-dimensional regions 
with boundary, |24j. This connection was first pointed out in Ref. [25J but without of- 
fering the details. Thus, it turns out that the renormalization group analysis of sigma 
models with embedded branes in their target space provide a natural field theoretic 
framework to address and interpret many important results that have been derived in 
the mathematics literature on the mean curvature flow. At the same time, new ideas can 
be brought in mathematics by implementing the perturbative and non-perturbative as- 
pects of quantum field theories with boundaries in the modern trends and studies. In this 
context, the mean curvature flow is tautonymous to the boundary renormalization group 
flow, to lowest order in perturbation theory, whereas the fixed points of the flow, which 



4 



correspond to conformally invariant boundary conditions, are the familiar D-branes in 
target space. Running solutions become applicable to the problem of tachyon condensa- 
tion in open string theory and to the Kondo effect of screening magnetic impurities by 
conduction electrons in metals. 

There has been considerable activity in recent years on boundary renormalization 
group flows and related quantum field theory problems, [26j [27J [28J, [29j [30j [31]. Most 
of the existing work concerns the construction of integrable quantum field theories of 
boundary interactions and their target space and world-sheet interpretation. These are 
based on mini-superspace models of the complete boundary flow, which is defined, in all 
generality, as dynamical system in the infinite dimensional space of all possible embedded 
configurations M in a given background Ai. The ambient space can be arbitrary, having 
its own renormalization group flow, but to simplify matters only models with boundary 
interactions on Ricci flat target spaces have been considered so far. Even the simplest 
case of boundary interactions in the quantum field theory of two free bosons, represented 
by embedded curves in R 2 , is quite rich and has not been fully explored yet. Specific 
proposals were made in this context for the exact form of the boundary quantum states 
of the so called semi-classical circular, paper-clip and hair-pin curves in R 2 , which take 
into account perturbative as well as instanton corrections, but their validity has only 
been tested in special limits. Boundary states associated to more general trajectories of 
the renormalization group flow do not seem tractable, to this day, due to the absence 
of systematic framework connecting the world-sheet with the target space description 
of deforming branes. The exact characterization of fixed points, satisfying appropriate 
Virasoro constraints, [32] , or (in some cases) extended conformal world-sheet symmetries, 
is also a difficult problem that has not been brought to the same level of understanding 
as for the bulk conformal field theories. Finally, the existence and construction of non- 
trivial infra-red fixed points, which take into proper account non-perturbative effects, as 
in the case of sigma models with 6 = n terms, have not been investigated in all generality 
(see, however, Ref. [26j-|31j for results in some special cases). Thus, it is only fair to say 
that the subject of boundary interactions and associated flows in quantum field theory is 
still at its infancy, in many respects, and any new insight cannot do less but contribute 
further to its development. 

Here, we will concentrate entirely on the target space description of running branes, 
as they arise in the semi-classical regime, examine general features of the mean curvature 
flow and obtain exact solutions. Although several of these solutions are known in the 
mathematics literature, there has been no proper mention or use in quantum field theory 
apart from some notable cases. Apart from bridging this gap, new solutions will also be 
constructed and studied in detail. Most results will be limited to two dimensions, where 
the mean curvature flow assumes its simplest - yet quite non-trivial - form for embedded 
or more generally immersed curves in R 2 . However, generalizations to curved spaces in 
two or higher dimensions will also be in focus, in which case the metric of the ambient 
space may also deform according to the Ricci flow. This combination of intrinsic and 
extrinsic curvature flows is quite natural from the physics point of view, as derived from 
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the generalized system of beta function equations for non-conformal sigma models with 
non-conformal boundary conditions. It should be contrasted to the form of the mean 
curvature flow appearing in the mathematics literature, where the metric of the ambient 
space is usually fixed once and for all. 

Further generalizations entail the inclusion of anti-symmetric tensor fields and Abelian 
gauge fields, which substitute the area functional of the mean curvature flow by the 
corresponding Dirac-Born-Infeld action, [33], [31], [35]. The presence of fluxes supports 
the existence of new solutions and alters the structure of the fixed points, as in the case 
of Ricci flow. In this context, the beta function of the gauge field also finds its proper 
place in mathematics in terms of the so called (Abelian) Yang-Mills flow. Thus, the 
sigma model approach to the closed and open sectors of string theory provide a unifying 
framework for studying all these different kind of geometric flows. In the general case 
one has to consider the coupled system of beta function equations for all massless modes 
of the string. The addition of fluxes will be addressed properly in future works. 

The presentation of the material is organized as follows. In section 2, the theory of 
Dirichlet sigma models is summarized and their boundary renormalization group equa- 
tion is identified with the mean curvature flow. In section 3, the flow is considered on 
the two-dimensional plane and various forms are derived in view of the applications. En- 
tropy functionals and their relevance to the structure of singularities are discussed. In 
section 4, several running solutions on the plane are introduced and studied in detail; 
they include special curves that evolve by translations, rotation or scaling, as well as 
other configurations that provide consistent mini-superspace truncations of the general 
evolution equations. Some of these solutions have already appeared in the physics litera- 
ture, but several others, like the Abresch-Langer curves (to be discussed later), have not 
yet found a proper place; they may serve as (p, q) models for boundary interactions in an 
appropriate setup. In section 5, a thorough analysis of the (in)stability modes associated 
to special solutions is performed and general results on the eigenvalues of the linearized 
operators are obtained in terms of supersymmetric quantum mechanics. Then, geometric 
transitions between different configurations are envisaged for various curves. In section 
6, generalizations of the mean curvature flow are considered on two-dimensional curved 
spaces. The simplest examples are provided by curves embedded in the Euclidean black 
hole background, with the familiar cigar shape, which is a Ricci soliton. Other examples 
include curve shortening problems on simple backgrounds that also deform by the Ricci 
flow, such as the sphere and its axially symmetric sausage-like variations. In section 7, 
we consider the mean curvature flow of surfaces in R 3 and then specialize to cylindrical 
branes as well as branes of revolution. Their evolution is reduced to an effective curve 
shortening problem on the plane. Finally, in section 8, we present the conclusions and 
list directions for further work. 

There are also three appendices included at the end for completeness. Appendix 
A summarizes the embedding equations of hypersurfaces in Riemannian geometry and 
provides the appropriate definitions. Appendix B compares the mean curvature flow to 
other systems of evolution equations for planar curves and draws connections to integrable 
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systems. Appendix C reviews the emergence of the mean curvature flow from magneto- 
hydrodynamics by dimensional reduction of the resistive diffusion of force-free magnetic 
fields. Several exact solutions discussed in the paper can be re-interpreted in this context 
and enjoy astrophysical applications. 

Parts of this paper can be considered as review of the main mathematical results on 
the mean curvature flow, but even in those cases there are supplementary details and 
alternative viewpoints that are of interest in physical applications. We hope that their 
systematic presentation will prove useful in many respects. The world-sheet description 
of various solutions, the role of non-perturbative effects in the characterization of the 
exact boundary states, as well as the addition of fluxes will be left to future publications. 

2 Dirichlet sigma models 

Consider a two-dimensional sigma model associated to maps of a two dimensional surface 
E into a general Riemannian manifold A4 with local coordinates and metric G^ U {X). 
It serves as classical model for string propagation in target space A4 of dimension m, 
whereas E is the two-dimensional world-sheet that is generally assumed to have boundary 
<9E. For practical purposes, E is taken to be a disc with connected boundary <9E = S 1 ; 
more complicated world-sheets, such as the annulus, may also be considered if higher 
loop string corrections are to be included in the study. Furthermore, (some of) the 
target space coordinates X 11 are taken to satisfy Dirichlet boundary conditions, i.e., 

X» | ra = fV) , (2.1) 

where y A are local coordinates in an n-dimensional submanifold Af of AA . These boundary 
conditions ensure that the variations of the sigma model fields along the world-sheet 
boundary, 5X^ \gs, are tangent to Af. According to this, there are m — n Dirichlet 
conditions imposed on the fields, in general, thus defining Af as classical brane embedded 
in AA. Then, the embedding equations (12 .ip follow the general theory of Riemannian 
geometry as outlined in appendix A. The extrinsic curvature of the classical branes Af, 
as well as the Ricci curvature of the ambient space AA, are arbitrary at this point, but they 
will be shortly constrained by quantum mechanical requirements if conformal invariance 
on E is to be maintained to lowest order in perturbation theory. 

2.1 Ricci and mean curvature flows 

With these explanations in mind, the corresponding Dirichlet sigma model is defined by 
the classical action 

S = -^— I d 2 z G^{X)5 ab d a X»d b X» + J—Idr G» v {X)V»{X)d n X» (2.2) 
Ana' js Zna' JdT, 

using conformally flat coordinates on E. Here, r denotes the parameter along the world- 
sheet boundary and d n is the derivative operator normal to it. The boundary contribution 
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to the usual two-dimensional action allows for the coupling of arbitrary vector fields in 
target space that are perpendicular to the submanifold Af. Thus, as it is customary, 
G>i,(X) is considered as generalized bulk "coupling constant" of the two-dimensional 
sigma model and V^(X) as the corresponding generalized boundary "coupling constants" . 
There are m — n independent vector fields of A4 that can couple to the normal 
derivatives of the fields X^ at <9£. Note the additional possibility to consider n fields 
coupled to the tangent derivatives of the coordinate fields y A in Af along the world-sheet 
boundary dH. They naturally form the components of a U(l) gauge field that lives on 
the brane; together with the anti-symmetric tensor field they may be used to provide 
flux generalizations of the present framework. Here, these additional fields are set to 
zero, thus only considering pure metric sigma models with embedded branes of arbitrary 
codimension. Throughout the paper, the signature of the world-sheet and of the target 
space will be assumed Euclidean. 

Two-dimensional sigma models, with or without Dirichlet branes, are perturbatively 
renormalizable quantum field theories. First, there is the bulk renormalization group 
flow of the target space metric that follows from the standard computation of the metric 
beta function, [2], 

P(G^) = Kfu, (2.3) 

to lowest order in perturbation theory. It is valid when all components of the curva- 
ture are small for otherwise higher order curvature terms arising from higher orders in 
perturbation theory become increasingly important. Such corrections will be excluded 
here, thus taking the lowest order result at face value. According to this analysis, sigma 
models are not scale invariant, in general, since their target space metric depends on the 
energy scale of the quantum field theory. In particular, identifying the logarithm of the 
world-sheet scale with the deformation variable t, one obtains the bulk renormalization 
group equation 

BC 

= -R, u + V^, + , (2.4) 

by also including the possibility to perform arbitrary reparametrizations along the flow 
generated by a vector field As such, it coincides with the general form of the 
(un)normalized Ricci flow for G^X; t). The fixed points of this flow, R^ u = V^+V^, 
are in accordance to the scale invariance of the world-sheet theory for all vector fields £ M . 
Note that the renormalization of the bulk metric is inert to the existence of embedded 
branes in target space. 

A special instance of these equations arises for gradient vector fields, £ M = — <9 M $, in 
which case $(Jf) assumes the role of the dilaton associated to anomalous transformation 
law of the target space coordinates of sigma models, <5 e X M = e<9 M $, under Weyl transfor- 
mations of the world-sheet metric, S t j a b = e^ a b- Recall, in this context, that the dilaton 
field $(Jf) enters into the bulk sigma model action as 

Shnlk = 4^L dh ( G ^( x h abd * xlld i> xu + dRhMX)) (2-5) 
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in order to ensure that the two-dimensional theory will remain renormalizable if the 
world-sheet metric 7 is not flat, having non- vanishing curvature R[y]- Then, the dilaton 
has its own beta function which together with the metric beta function yield the modified 
system of renormalization group equations of the bulk theory, 

J^G> = -f3{G»v) = -R»» - 2V M V,$ , (2.6) 

1* = _/?($) = -(^$)(^$) + -V 2 $ + -^-j- . (2.7) 

The last term above accounts for the central charge of the model and it can be arranged 
so that it cancels in critical string theory. Here, the central terms will be kept, since the 
dimension m is arbitrary. 

In any case, these beta functions satisfy m differential identities derived from the 
non-renormalization condition of the trace of the energy-momentum tensor of the sigma 
model, [36J, 

<y?($) = (d»<S>)(3(G^) + -V (>(GV) - l -G, v G xp P{Gx P )) , (2.8) 

which is valid for all t to lowest order in a'. Weyl invariance of the two-dimensional 
theory is only achieved at the fixed points of the modified (G^, $) flow; it should be 
compared to the weaker condition of scale invariance that was considered earlier. At the 
fixed points of the background metric flow, j3{G fa ,) = 0, equation (12.81) implies that /3($) 
is constant on M. and can be set equal to zero without loss of generality. Thus, Weyl 
invariance implies j3(G^ u ) = = /?($) simultaneously. Finally, note that the dilaton flow 
is also inert to the existence of branes in target space, as for the metric. 

For sigma models with branes in their target space one also has to consider the 
dependence of the embedding equations (12. ip on the energy scale of the two-dimensional 
quantum field theory. This calculation was first performed in all generality in Ref. |24j, 
where it was found that the deformations of Dirichlet branes, as described by the one- 
parameter family of functions f^(y A ; t) that may depend on the logarithmic world-sheet 
length scale t, are driven by their extrinsic curvature vector to lowest order in perturbation 
theory. The analysis is performed by first considering the variation of the classical action 
( 12.21) leading to the following set of compatible boundary conditions, 

f^G^dnX" = , V = , (2.9) 

with respect to the embedding equations of N into M.. If one expands around such a 
reference configuration, denoted by X M , will decouple completely from the compu- 
tation. However, if one considers quantum corrections, there will be counter-terms for 
V 1 that can be computed by regulating divergences of the relevant graphs with a short 
distance cutoff e. Starting with the reference configuration X M , with X M = f^(y A ) on 
<9E and introducing normal coordinates in target space and on the brane, as usual, the 
counter-term assumes the form 

AS = ~ h L dT G " v (? ABK *b<) d nX u (log e) , (2.10) 



9 



to lowest order in perturbation theory, and changes V th accordingly. 

Thus, to this order, the associated beta function for the boundary coupling is 

/8(V") = -g AB K° AB h» . (2.11) 

Here, and above, the right-hand side involves the trace of the second fundamental form 
K\ B of the brane with respect to the induced metric g on it, is a complete basis of 
unit normal vectors to the brane in Ai and a labels the transverse directions. In turn, 
one arrives at the following boundary renormalization group equation for the embedding 
conditions (|2.ip . 

-j- = H^-e, (2.12) 

using the notation of appendix A for the mean curvature vector normal to the brane. 
Here, we have also included the freedom to perform arbitrary reparametrizations along 
the flow generated by a vector field it is the same vector field that enters into the 
Ricci flow (12. 41) . The resulting equation coincides with the general form of the (unnor- 
malized) mean curvature flow studied in mathematics. We will take this equation at face 
value and suppress all possible higher curvature terms3 which may arise at higher orders 
in perturbation theory. Within this approximation, and from now on, the two terms 
"boundary renormalization group flow" and "mean curvature flow" will be used without 
distinction. 

The dilaton $ is introduced by generalizing the Dirichlet sigma model (12.21) to curved 
world-sheets, so that the two-dimensional action consists of the bulk term (12. 5p plus 
boundary contributions, 

S = S bulk + -^/ dr(G flI/ (X)V^(X)d n X v + a , K^(X)), (2.13) 

thus also taking into proper account the coupling of the dilaton to the extrinsic curvature 
K of the boundary. The quantum theory is renormalizable and the boundary flow is 
provided by equation (I2.12p above with £ M = — <9 M $. However, it should be noted at this 
point that if the ambient space exibits isometries generated by a Killing vector field 
the choice = — <9 M $ + will affect the mean curvature flow but not the Ricci flow. 
Recall that d^k u ^ vanishes identically and so (3{G^ V ) does not change. As for the dilaton, 
it can be consistently taken to satisfy the relation /c M <9 M $ = and also does not 

change. Thus, apart from the standard bulk flows, we obtain the following boundary 
flow 



ot 



H a h% + - fc" . (2.14) 



Weyl invariance of the quantum Dirichlet sigma model leads to fixed points of the 
combined buck and boundary renormalization group equations satisfying the general 



J Such terms have been computed systematically in some cases in Ref. [37] and [35], in analogy to 
higher curvature correction terms computed for the Ricci flow, [2], on general grounds; we thank Arkady 
Tseytlin for bringing some of those references to our attention. Further results can also be found in the 
more recent work Ref. 139 . 
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relations 



(2.15) 



supplemented by the vanishing condition for the dilaton beta function, when $ is non- 
trivial, 



Fixed points of this type will be discussed later in section 6. 

According to all this, Dirichlet sigma models provide a natural framework to real- 
ize and unite both Ricci and mean curvature flows, since one has bulk and boundary 
renormalization group equations defined with respect to the same deformation variable 
t. Simple fixed points correspond to backgrounds with Ricci flat metrics in which there 
are embedded branes as minimal submanifolds (of arbitrary codimension), so that their 
extrinsic curvature vanishes; these are the familiar D-branes. More general fixed points 
also arise in the presence of non-trivial dilaton field. In either case, the corresponding 
solutions are associated to two-dimensional conformal field theories defined on a disc 
with conformally invariant boundary conditions. Away from the fixed points one has 
deformations of branes in deforming metric backgrounds, in general, but there is also the 
simpler possibility to consider Dirichlet branes with non-conformal boundary conditions 
deforming in backgrounds with fixed metric satisfying bulk conformal invariance. The 
simplest example of this kind arises in the two-dimensional quantum field theory of sev- 
eral free bosons, in which case M. = R m , and impose non-conformal boundary conditions 
on <9E so that the branes will not be embedded as minimal submanifolds. Even in such 
simple cases there can be many interesting possibilities, as will be seen later; also the 
systematic construction of the corresponding boundary states in quantum field theory is 
far from being complete, up to this day. 

2.2 Gradient flow description 

It is well known fact that the mean curvature flow of branes can be formulated as gradient 
flow of their volume functional, 



given in terms of the determinant of the induced metric Qab- Indeed, simple calculation 
shows that the first variation of the volume with respect to the embedding variables / M 
yields 




26 — m 
6a' 



(2.16) 




(2.17) 




(2.18) 



In the presence of dilaton one has to consider the effective volume functional 




(2.19) 



and derive the generalized mean curvature flow as gradient flow 




(2.20) 



11 



with 

^ - T^m < 2 - 21 > 

which is positive definite. V[/, $] can be extended to the full Dirac-Born-Infeld action 
in the presence of fluxes, |24j . 

When an evolution equation arises as gradient flow of the general form 

it is natural to investigate its dissipative character and try to associate with it monotonic 
functionals in time. S[ip] itself evolves in time as 

d qUn] OS dip 1 GlJ dS dS 

Jt s[ip] = W^ = ~ G WW { ] 

and, therefore, it decreases along the flow when Q is positive definite. The mean curvature 
flow is an example of this kind, and, naturally, the branes deform by lowering their total 
volume towards minimal submanifolds. Of course, there can be other functionals which 
are also decreasing monotonically in time and serve as entropy of the deforming data. 
Examples of this will be encountered in section 3 for the mean curvature flow defined 
in flat ambient spaces. In general, there is no straightforward procedure to construct 
entropy functionals for gradient flows, which may exist irrespective of the positivity of 
Q; see, however, [30], for some recent general results in this direction. 

Similar considerations can be applied to the Ricci flow for comparison. The Ricci flow 
arises as gradient flow from the Einstein-Hilbert action 

S E [G] = f d m X Vd^G R[G] (2.24) 

J M 

i.e., ^ 

t) = W^3fJ = ■ ( 2 - 25 ) 

In this case, the appropriate matrix Q is provided by the DeWitt metric in superspace 
consisting of all target space metrics on M., 

G^x5G^5G kX = \J m d m X v/detG (sG^SG^ - \{5G%){5G\)^ . (2.26) 

In the presence of dilaton, the modified Ricci flow is also described as gradient flow 
using the Einstein-Hilbert-dilaton action 

S E [G, $] = / d m X VdrtG e" 2 * (r[G) + 4(<9,$)(d /i $) + 2 26 ~™ ) . (2.27) 

More precisely, setting ip 1 = (G^, $), one finds that the beta functions of the metric and 
dilaton fields take the form 

V = _ g ijSSM ; SS\<P\ = _g u ^L (2 . 28) 
dt 5ip J 'Sip 1 dt 
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with 



/ 4:Gfj,\G u p G^y \ 

V G Xp \{m-2)J 



(2.29) 



and 



/ ~(G^G up - \G^G x f) G» v \ 
v G Ap -4 / 



(2.30) 



The latter expression generalizes the DeWitt metric to the metric-dilaton system so 
that ( 12.26P is replaced by 

Qu&<ft8<p J = - f d m X VdetG e" 2 * (sG^GT - -{5G\ - 45$) 2 ) . (2.31) 
4: Jm V 2 r J 

Remarkably, the Einstein-Hilbert-dilaton action is also the total time derivative of the 
effective volume of Ai, 

S E [G, $1 = 4 (7 d m X e" 2 * V&itG) . (2.32) 
at \Jm J 

Note, however, that the DeWitt metric is not always positive definite for it is well 
known that the Weyl mode fluctuations of the metric on Ai have negative norm and 
they are naturally associated to "time-like" directions in superspace. By the same token, 
the extended DeWitt metric on the metric-dilaton superspace also exhibits "time-like" 
directions which now arise from a combination of the dilaton and the Weyl mode of the 
target space metric. Thus, unlike the case of mean curvature flow, Se does not vary 
monotonically with time. However, there is a closely related functional introduced by 
Perelman, which serves as entropy for the Ricci flow on compact Riemannian manifolds. 
It can be thought as being inspired by string theory constructions combined with a special 
choice of reparametrizations along the flow, but the full details are beyond the scope of 
the present work. We only mention here that 

\[g]:= mm W S E {G,<5>) with f d m X e" 2 * VdrtG = 1 (2.33) 

provides a monotonically increasing functional along the Ricci flow, [10] , by removing the 
unwanted "time-like" directions of the DeWitt metric. Then, X[g] is interpreted as the 
lowest eigen- value of the operator —V 2 + -R/4, which is defined in terms of the metric at 
any given moment t, whereas the constraint on the effective volume of Ai provides the 
normalization of the corresponding eigen-function exp(— $). In effect, X[g] is determined 
by applying the variational method of elementary quantum mechanics on Ai. A simpler 
version of this construction appeared first in the physics literature, [H], and enters into 
the definition of the (monotonically decreasing) effective central charge along the Ricci 
flow; see also Ref. [12] for more details and further generalizations, as well as Ref. [13] for 
its extension to higher orders in a' and in connection with Zamolodchikov's c-function, 
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[14"] . Other important entropy functionals have also appeared in the literature, [10J, but 
their physical interpretation is not yet as clear. 

All properties and entropy functionals of the bulk flows are independent of the exis- 
tence of embedded branes. On the other hand, the analysis of boundary flows depends 
crucially on the metric of the ambient space. Important entropy functionals for the mean 
curvature flow will be discussed later. It should be noted, nevertheless, that most results 
in mathematics are concerned with the mean curvature flow in flat Euclidean space or 
in curved Riemannian spaces with fixed metric, apart from some notable exceptions^]. 
The physical origin of these flows suggests that they should be studied together in all 
generality and new entropy functionals should be found. 



3 Mean curvature flow on the plane 

The simplest framework for studying the boundary renormalization group flow of Dirich- 
let sigma models is provided by the two-dimensional quantum field theory of two free 
fields whose values at the boundary of the world-sheet are restricted to lie on a given 
curve. In this case the target space is R 2 and trivially satisfies conformal invariance for 
the bulk metric beta function. Thus, the only interesting thing to consider are bound- 
ary effects, which in general are associated to deformations of the Dirichlet curve due 
to renormalization on the world-sheet. Since there is no dilaton or any other additional 
fields in this model, the boundary renormalization group flow is identical to the mean 
curvature flow of embedded curves in the plane, 

dr -> 

s = ™-f. (3-D 

where f is their position vector and £ includes the effect of reparametrizations along the 
flow; immersed curves can also be considered by allowing for self-intersections. 

Although £ will be left arbitrary in the mathematical presentation below, conformal 
invariance of the free field theory in the bulk requires that it can only be a Killing 
vector field on the plane so that the target space metric remains at its trivial fixed point; 
otherwise, one has to reabsorb it into the time evolution of the curve and eliminate it 
all together. Turning on a general £ may lead to mathematical simplifications of the 
curve deformations; of course, the tangential part of the deformations can always be 
removed by appropriate diffeomorphisms. In either case, the mathematical structure of 
the equation is the same although the physical interpretation of its solutions differs. We 
will always insist on having conformal invariance for the bulk space theory and only allow 
for non-conformal boundary conditions. 

2 Some aspects of the combined system of Ricci and mean curvature flows have been considered by 
Hamilton, [45) : we thank Klaus Ecker for bringing this to our attention. 
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3.1 Basic general elements 



Let us first consider various forms of the mean curvature flow for embedded curves (open 
or closed) in the plane, which are convenient for later use and also help to set up the 
notation. The points of R 2 are parametrized by the position vector r with Cartesian 
coordinates (x,y) and any given curve will correspond to an orbit f(s) = (x(s),y(s)) 
with respect to an afhne parameter s. Alternatively, one may think of a curve as the 
graph of a function y = <f(x) when x is identified with s. Such curves are not stationary 
but they evolve according to the mean curvature flow with respect to the deformation 
time t so that the corresponding trajectories are parametrized in Cartesian coordinates 
as (x(s,t),y(s,t)) or in equivalent graph form as y = ip(x(t),t). 

The tangent vector at each point of the curve is drj ds and therefore the unit normal 
vector inward to the curve is 

ft 1 ( J_V 9x\ _ (-y'Qr), 1) 

^{dx/dsf + (dy/ds) 2 V ds' ds) j\ + ^{x) ' 

Furthermore, since the induced metric (line element) on the curve is 

.*=((£Y+(£YW. (3.3) 



where I is the arc-length (or proper length) on the curve, it follows by definition of the 
mean curvature H that 

H= 1 ( d 2 ydx d 2 xdy\ _ <p"(x) 

The arc-length of the curve can be used to cast the mean curvature flow in the form 
df/dt = d 2 f/dl 2 , which resembles the heat equation, albeit is non-linear, but this is 
not very practical for finding explicit solutions. Instead, the mean curvature flow in R 2 
assumes the following convenient form, also taking into account arbitrary reparametriza- 
tions generated by a vector field £ along it, 

dt 5 (i + v? /2 ) 2 ' dt~ * (i + ^ 2 ) 2 ' 1 j 

Note at this point that since y(t) = (p(x(t),t) we have 

9y dl P . / / \ dx 

m = m +v{x) m> (3 - 6) 

which in turn implies the following simple form of the mean curvature flow of graphs 
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The fixed points are characterized by the second order equation 

(arctany/(x))' = ? - <p'C (3.8) 

that also includes the effect of arbitrary reparametrizations in their classification. 

Another convenient form of the mean curvature flow in R 2 follows by considering the 
mean curvature H as function of the slope of the curve, 

/3 = arctan<//(x) , (3.9) 

which is the angle formed by the tangent at each point of the curve with the x-axis. In 
terms of this variable, the unit (inward) normal vector is n = (—sin/3, cos/3), whereas the 
unit tanget vector is t = dr/dl = (cos/3, sin/3) at each point. Considering the projection 
of the position vector r onto the unit normal, S(f3) = —r(/3) ■ h, it follows that the 
Cartesian coordinates of the curve can be expressed as functions of the slope, 

x (/3) = S'(/3)cos/3 + 5(/3)sin/3 , y((3) = S'(/3)sin/3 - 5(/3)cos/3 , (3.10) 

where prime denotes derivative with respect to (3. Furthermore, we have the identity 

S{/3) + S"((3) = x'(/3)cos/3 + ?/(/3)sin/3 = — • t = (3.11) 

since dl/d/3 = l/H(/3). Then, upon differentiation of equations (I3.10p . it turns out 
that x'{0) = cos/3/if(/3) and y'(/3) = sin/3/if(/3). As a result, the curves are fully 
determined, up to translations, by specifying the mean curvature H as function of the 
slope /3, according to the relations 

When the curves deform by the mean curvature flow, with £ = 0, the evolution for 
S((3,t) satisfies the simple relation dS/dt = —H that follows from its definition. Then, 
employing the identity (13.111) . one easily finds that H({3,t) satisfies the parabolic partial 
differential equation 

9H tf 2 S + ^- (3-13) 



dt " a/3 2 



This form will be particularly useful for understanding the characteristic features of some 
special solutions listed in section 4. When reparametrizations generated by £ are also 
included along the flow, one finds dS/ dt = —H + £ • n and the differential equation for 
the mean curvature generalizes to 

f = H w +H *- H2 (w (i - h)+i - h )- (3 ' 14) 

For locally convex closed curves with winding number n, there is a periodicity con- 
dition on both coordinates, x(/3 + 2nn) = x(j3) and y(/3 + 27m) = y(/3), which implies 
that 

2-ku piP 

m d P = . (3.15) 
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The extrinsic curvature of such curves satisfies the periodic condition H(/3-\-27in) = H(j3), 
but there can be cases of closed curves, as will be seen later, where H(/3) has smaller 
period. Also note that solutions with periodic extrinsic curvature do not necessarily yield 
closed curves no matter how many times they are iterated. An elementary example of this 
kind corresponds to the choice H(0) = 1 + cos/3; it yields the curve x{(3) — (3 — tan(/9/2) 
and y(fi) = -log(cos 2 (/3/2)) so that x(f3 + 2vr) = x(f3) + 2vr and y(J3 + 2vr) = y(j3). 
Likewise, for H{(3) = 1 + sin/3 one has x(/3 + 2tt) = x(f3) and y({3 + 2n) = y(/3) + 2ir. 

Finally, we also include for completeness the form of the mean curvature flow on the 
plane using polar coordinates x = rcosO and y = rsin#. In this case an arbitrary curve 
on the plane can be thought as graph r = p{6) that evolves in time according to 

d P _ idp_ p i + 2p> 2 -pp" 
dt pde p{p 2 + p' 2 ) ' 

Here, prime denotes the derivative with respect to 9 and the evolution of p(8(t),t) is 
computed using dr(t)/dt = dp(6(t),t)/dt + p'(6)d6(t)/dt. Arbitrary reparametrizations 
along the flow can also be included, if needed, in the system of polar coordinates. 



3.2 Entropy functionals, curvature bounds and singularities 

The mean curvature flow df/dt = Hh tends to deform curves in the direction of their 
inward normal vector, as if there were tension forces depending on the magnitude of H 
at each point. As a result, open lines tend to become straight, whereas closed curves 
tend to become round circles as depicted in Fig.l. 



(a) (b) 

Figure 1: Evolving open and closed curves on the plane 

The area A surrounded by closed curves 7 also tends to decrease at constant pace as can 
be easily seen by computing 

f = />* = -/J^-/„> = - 2 - w 

Thus, closed curves have the tendancy to shrink, as they become rounder and rounder, 
[IB]) [S]) until they fully collapse to a point at some time T, which for all practical 
purposes can be taken to be zero. 
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Based on this observation one may consider rescaling the coordinates and redefining 
time as 

Hs,t) = M) , t = --log(-2t) , (3.18) 

so that the evolution takes the equivalent form 

JLr(s, t) = #(s, t)n + r(s, t) (3.19) 

in terms of the mean curvature H of the rescaled curve in R 2 ; it so happens that 
H = y/—2tH. The variant (I3.19P is called normalized mean curvature flow, since, by 
construction, it preserves the area surrounded by the rescaled closed curves with respect 
to the new time variable t. As — oo < t < 0, we see that — oo < t < oo and so the 
normalized solutions exist for all time. The normalized flow can be alternatively viewed 
as special instance of the unnormalized flow (13.11) when reparametrizations are performed 
along it with £ = — r and the tilde is dropped for comparison. This will be quite use- 
ful later for understanding the structure of scaling solutions and the characterization of 
singularities that may form by the flow. 

There are entropy functionals associated to the mean curvature flow. First, let us 
consider the backward heat kernel, defined on R 2 for all t < 0, 

/C(r, t) = ^— exp ( -^j* ) > (3- 2 °) 

and integrate it over the curve by the induced arc-length that also varies with time. It 
follows that this is a monotonically decreasing functional, due to Huisken, [15J , since 



d 
dt 



j K,{r,t)dl = - J K,(r,t) (h + dl < (3.21) 



by the unnormalized flow df/dt = Hn which is applied to closed curves C. Since S = 
—r ■ h, any closed curve that satisfies the special relation 

(-2t)Hh + r = (3.22) 

keeps the entropy functional invariant. Configurations of this kind are self-similar solu- 
tions with factorized time dependence. The simplest example is provided by a uniformly 
shrinking round circle whose radius varies as \J—2t and H as 1/ \/—2t. 

Likewise, for the normalized mean curvature flow, we consider the Gaussian weight 
function _ 2 

K(r) = -Lexp ( - T — j (3.23) 



. 2vr V 2 , 

that depends implicitly upon t, and integrate it over the rescaled curve C by the corre- 
sponding induced arc- length I. It follows again that this is a monotonically decreasing 
functional with respect to the rescaled time, [iS| . since 

4 / £{r)dl = - { iC{r){H - S) 2 dl< . (3.24) 
at Jc Jc 



18 



In this case the extrema of the entropy functional satisfy the normalized self-similar 
condition 

Hh + r = 0, (3.25) 

which is attained as t — > oo. The round circle is a fixed point of the normalized flow. 
Other non-trivial fixed points also exist, and they are classified by the so called Abresch- 
Langer closed curves, as will be seen in the next section. 

Next, we discuss certain bounds on H that are important for the classification of 
singularities formed by mean curvature flow. There is always a lower bound for the blow- 
up rate of the curvature which is derived by applying the maximum principle. Indeed, 
specializing equation (13.131) to the maximal value H max _(t) attained at each instant of 
time, we obtain the inequality 

Q 3 3 

^^max — H max ^ H max + -ff max < -f^max • (3.26) 

Closed convex curves have H max > and they develop curvature singularities at some 
finite time, say T = 0. More generally, it follows 

H max {t) > , (3.27) 

by integrating the inequality from t to 0, thus establishing a universal lower bound for all 
t < 0. The uniformly shrinking round circle saturates this curvature bound at all time. 

On the other hand, it is not at all guaranteed that there is an analogous upper bound 
for H max (t) based on general grounds. Actually, the singularities of mean curvature 
flow are divided in two general categories. Their characterization is based on bounds of 
\H\ max (t), which is taken in absolute value in general. Type I singularities arise when an 
upper bound of the following form also exists, 

> \H\ max (t) > -jL= , (3.28) 



with appropriately chosen constant C < oo. This is equivalently stated as 

C > \H\ max (t) > 1 (3.29) 

using the rescaled curvature H. In all other cases the singularities that are formed are 
called type II. All closed embedded curves in R 2 will eventually form type I singularities. 
Even if the curve is not convex at a given time it will become convex at later times 
and follow the evolution towards the singularity by becoming rounder and smaller, i.e., 
asymptote the uniformly contracting circle, [16], [H]. However, one can also imagine 
deformations of immersed curves, with self-intersections, whose curvature blows up at 
faster rate as they begin to develop cusps and yield type II singularities. 

All closed planar curves satisfying the curvature bound (13. 28ft tend to self-similar 
solutions (13.221) in the vicinity of the singularity, [H] (but see also Ref. [19] and [50]). 
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The proof relies on the monotonic behavior of Huisken's functional whose extrema are 
the self-similar solutions. Thus, the classification of type I singularities reduces to the 
classification of self-similar solutions, which are completely known on the plane. The 
circle is the only embedded curve of this kind whereas the other solutions are special 
self-intersecting convex curves. In higher dimensions one may consider hypersurfaces 
that evolve by their mean curvature vector in flat space and generalize the constructions 
and results mentioned above. However, there is no systematic classification of the self- 
similar hypersurfaces that extremize the corresponding Huisken functional unless the 
hypersurface is compact with positive mean curvature. As a result, the general structure 
of type I singularities is less understood in higher dimensions. 

Type II singularities arise when the curvature blows faster than 1 / \/—2t as t — > 0~. In 
this case there is a sequence of times t n — ► CT such that the curve obtained by appropriate 
magnification at each instance t n , so that its maximal curvature becomes 1, will converge 
to a translating solution, [19] (but see also Ref. [50]). The latter is a very special solution 
of the mean curvature flow that will be discussed later together with the scaling (self- 
similar) solutions on the plane. A typical example of this kind of singular behavior arises 
from the evolution of a cardioid. It is a convex closed curve with winding number 2 that 
self-intersects once and consists of two loops, the inner and outer, touching each other. It 
can be intuitively seen that the inner loop will contract faster and form a cusp before the 
outer loop has a chance to shrink to zero size. Then, as one zooms closer and closer to 
the diminishing inner loop, as it begins to form a cusp, the shape of a translating solution 
will emerge according to the general statement above. The basic idea is illustrated in 
Fig. 2 below focusing on the rescaled region of maximal curvature. 



Similar results apply to the structure of type II singularities of evolving compact 
hypersurfaces of positive mean curvature in flat space. 

For convex planar curves 7 there is an additional entropy functional which is defined 





as follows, 




logH . 



(3.30) 
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It can be easily seen following, for instance, [18] that 



d 2 S „(fd£\ 2 dS 



using the normalized variant of the flow that exists for all t < oo. Note that if dS/dt 
were positive at some time, it would blow up at later times because d 2 £/dt 2 would also 
be positive. However, this is impossible for £ exists for all values t. Thus, we conclude 
that 

§<0 (3.32, 

for all t < oo. The round circles extremize £. 

This functional is the mean curvature analogue of a similar expression 

H = J cfxVdetG R hgR (3.33) 

introduced for the Ricci flow on two dimensional surfaces, [51]. Using the normalized 
Ricci flow on compact surfaces with R > 0, which are analogous to convex planar curves, 
it follows that Ti, decreases monotonically. Furthermore, the curvature remains bounded 
for all time, as 

C > R > c> , (3.34) 

with appropriately chosen constants c and C < oo. These bounds are analogous to (13.291) 
for the normalized mean curvature flow. Note that the round spheres extremize 7i. 

Finally, another important result on the subject states that the number of self- 
intersections of immersed planar curves can not increase by the mean curvature flow 
in the forward time direction, 



It will be interesting to explore new methods for integrating the mean curvature flow 
on the plane, at least formally, and devise a Lax pair formulation for it, if it is at all 
appropriate, by developing analogies with the algebraic treatment of various intrinsic 
curvature flows on two dimensional surfaces, [53J. Until then, we can only rely on the 
general mathematical results concerning the qualitative behavior of the flow, as outlined 
above, and the construction of various explicit solutions that will follow next. Note at this 
end that there are other type of evolution equations for planar curves leading to known 
integrable systems, as explained in appendix B that is only included for comparison. 



4 Special solutions on the plane 

Several exact solutions of the mean curvature flow on R 2 are listed and comments are 
made relating their appearance in the physics and mathematics literature. Many details 
will be filled in for completeness and a number of new results will also be derived. We will 
speak in Euclidean terms calling a planar curve Dl-brane, as opposed to the Lorentzian 
version of D0-branes having one- dimensional world-volume. 
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Although boundary interactions in quantum field theory provide the main framework 
for our work, as in Ref. [2E]-[2I] where a few explicit solutions have been constructed, it 
should be noted that various running solutions were also constructed a long time ago in 
different physical context. They first appeared in the original work on the motion of grain 
boundaries in an annealing piece of metal, [H], and later in the magneto-hydrodynamic 
theory of resistive diffusion of force- free magnetic fields, |54| . which is summarized in 
appendix C. Frequent references will be given to them at appropriate places in the text. 
In mathematics, some of these solutions are briefly discussed in the textbook [XSj and 
references therein. However, not all of them have yet found their exact place in quantum 
field theory. 

4.1 Trivial fixed points 

The fixed points of the mean curvature flow (13.71) are simply described by dtp/dt = and 
(f"(x) = when no reparametrizations are taken into account (£ = 0). Obviously, these 
are time invariant straight lines of the general form 



whose extrinsic curvature vanishes identically. They represent Dl branes on the plane, 
which are compatible with conformal invariance. DO branes also arise as points with 
fixed position in the (x, y) plane; they can be thought as the end-point of shrinking 
closed planar curves. 

4.2 Translating solution 

The simplest static solution of equation (13. 7ft . modulo reparametrizations, corresponds 
to the choice of a translational Killing vector field along the y-direction, £ = d/dy with 
components £ x = and £ y = 1, up to a constant factor v. Then, it follows that the shape 
of the curve is given by the graph of the function 



where a and b are integration constants. Setting them equal to zero amounts to placing 
the tip of the curve at the origin of the coordinates (x, y), in which case case it asymptotes 
the lines x = ±ir/2v. 

The resulting curve is called grim-reaper in the mathematics literature or hair-pin in 
the physics literature, [27], [28], where it was encountered before. The same configuration 
was also found in the earlier works [14J and [51] (but see also Ref. [20] among others). It 
is a translating solution along the ^/-direction for it can be alternatively viewed as moving 
linearly in time along the ^/-direction with constant spead v (in appropriate units), i.e., 
dtp/dt = v, so that 






(4.2) 



v 




(4.3) 



v 
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The graph of this configuration is given in Fig. 3 below. 



y 



Figure 3: The hair-pin (or grim-reaper) curve on the plane 

A hair-pin facing in the opposite direction is obtained by setting y —>■ —y, in which 
case the corresponding time dependent configuration translates linearly in time towards 
the negative ^/-direction. Thus, the sign of v selects one of the two possible cases: the 
hair-pin or the anti-hair-pin. The slope of these curves is simply given by (3 — vx, in 
which case their mean curvature is 

H((3) = vcos(3 . (4.4) 

The solution satisfies the special condition H"(f3) + H(f3) = that follows from equations 
(I3.13P or (13.141) depending on how one views the evolution. Likewise, a translating 
solution along the x-axis, in either direction, arises by reversing the role of the x and y 
coordinates. 

In all cases, the hair-pin solutions represent Dl branes supported by a linear dilaton 
$ in the coordinate of the translating direction. For example, a hair-pin that moves 
with constant velocity v along the y-axis has £ M = — with $ = —vy. In this 
respect, the translating solutions are examples of mean curvature solitons that do not 
affect the conformal field theory on the plane. Also, according to the variational method 
of section 2.2, in the presence of dilaton, such configurations represent geodesies on the 
plane equipped with the metric 

ds 2 = e 2vy (dx 2 + dy 2 ) . (4.5) 

This is a Ricci flat metric as it relates to the Euclidean frame in polar coordinates by the 
change of variables 

vr = e vy , 9 = vx , (4.6) 

so that ds 2 = dr 2 + r 2 d9 2 with 1 < vr < oo and -n/2 < 9 < n/2 for v > 0. This 
maps to a domain on the right half-plane plane, which is exterior to the disc vr < 1, and 
the hair-pin corresponds to the vertical straight line vrcos9 = 1 that is tangent to its 
boundary. 

The hair-pin is a gradient soliton, which should be thought as the mean curvature 
analogue of the well known Ricci flow soliton associated to the intrinsic geometry of a two- 
dimensional cigar, [5Tj ; the latter has the interpretation of a two-dimensional Euclidean 
black hole in conformal field theory, [55]. In appropriate context, the hair-pin serves as 
model for studying tachyon condensation in open string theory (see, for instance, [56] 
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and references therein). When v — the configuration becomes straight line, y = 0, as 
viewed from the origin of coordinates; alternatively, when this limit is considered from 
the view-point of an asymptotic "observer", situated at y = oo, the hair-pin looks like a 
semi-circle with infinite radius. 

Finally, we note the important property of this solution to exhibit just one point of 
maximal curvature situated at its tip. This is not accidental but consequence of a general 
theorem stating that any strictly convex solution of the mean curvature flow that exists 
for all time — oo < t < oo and the mean curvature becomes maximum at only one point, 
must necessarily be a translating soliton, [57] . This is an important ingredient that goes 
into the study of type II singularities, which look like a hair-pin following a sequence of 
appropriate magnifications that keep the maximal curvature normalized to a fixed value, 
v — 1, through out the evolution. 



4.3 Rotating solution 

Using the rotational Killing vector field £ = d/d6 = —yd/dx + xd/dy on the plane, up 
to an overall constant factor u, other static solutions of (13.71) follow by integrating the 
differential equation 

(arctan^'(x))' = uj (x + y(p'(x)) . (4.7) 
Then, for a graph y = (p(x), one obtains 

arctan<//(a;) = ~^{ x2 + y 2 ) + c = + c (4-8) 

in terms of the polar coordinate r. Setting the integration constant c equal to zero 
amounts to placing the curve at the origin of the coordinate system so that it starts 
tangentially to the x-axis. Further integration of the equation can not be performed in 
closed form, but this is no problem for drawing the shape of the resulting curve. Its slope 
at each point, as it reads from equation ( 14. 81) . is proportional to the distance-squared 
from the center, and, therefore, it follows the shape of an unbounded spiral as depicted 
in Fig. 4. 



(a) (b) 
Figure 4: Yin- Yang curve on the plane 

Thus, the integral curve of the differential equation H4.8[) is the unique rotating solu- 
tion on the plane, which is another example of a mean curvature soliton also known as 
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Yin- Yang curve in the mathematics literature; in physics it has appeared much earlier in 
Ref. [2]. Half of it is plotted in Fig.4a, whereas the complete curve appears in Fig. 4b 
by also plotting its symmetrical part about the origin, i.e., its point of inflection. The 
turns of the spiral are separated by approximately 2n/uor for large r. An alternative 
(dynamical) description of the solution is provided by a revolving spiral with constant 
angular velocity u about the origin, with simple time dependence 

axct(m<p'(x(t), t) = ^{x 2 + y 2 ) + c = ^r 2 + tot , (4.9) 

without making use of the compensating Killing vector field that stabilizes its rotation. 
If polar coordinates were used on the plane, the rotating solution (14.91) would be the 
graph of a linearly evolving function r = p{9 — ujt); in this case, equation (I3.16P simplifies 
a lot since d/dt = —ud/d6 when acting on p. In this frame, the mean curvature of the 
solution specializes to H(j3,t) = H{(3 — ui) and satisfies the non-liner equation in (3 

H 2 ((3) (H"(J3) + H{p)) + uH'ffl = , (4.10) 

as follows from equations (13.131) or (13.141) depending on the view-point. It is a non-trivial 
task to find the solution in closed form. 

The rotating solution represents an evolving infinitely long Dl brane with one point 
held fixed at the origin. Its existence, as static solution of equation (13. 7p . is attributed 
to the Killing vector field £* = cu(—yd/dx + xd/dy) having V x £ = 2uj ^ 0. As such, 
£ can not be derived from a potential as gradient field, and, hence, there is no dilaton 
field associated to the static form of the solution that could account for its boundary 
conformal field theory description. Thus, unlike the grim-reaper, the Yin- Yang curve is 
not a gradient soliton. It corresponds to a boundary quantum state that does not satisfy 
conformal invariance, but it runs, via rotation, with the world-sheet energy scale. The 
sign of uj distinguishes the two different modes of rotation around the clock. 



4.4 Scaling solutions 

Another type of static solutions of equation (13.71) follow by considering the vector field 
£ = rd/dr = xd/dx + yd/dy, up to a constant factor c, which generates dilations 
of the plane. These are called scaling (or homothetic) solutions since they evolve by 
overall scaling when the time dependence is reinstated at the expense of suppressing 
the corresponding reparametrizations along the flow. As such, they satisfy the defining 
relation 

(axctanip' (x)) = c(<p(x) — x<p\x)) (4-H) 

for y = (f(x). Note, however, that the generator £ of dilations is not a Killing vector field 
on the plane, although £ M = — with 

®(x,y) = ~(x 2 + y 2 ) . (4.12) 
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If conformal invariance of the quantum field theory of the plane is to be maintained, 
the homothetic solutions will unavoidably arise as time dependent curves with factorized 
t-dependence so that 

y = ( p(x(t),t) = VM<p(-j=,l\ . (4.13) 



The sign of c determines the basic features of time evolution. Note that in all cases 
ct should be strictly non-negative. Thus, for c < 0, the scaling solutions are shrinking 
as t runs from — oo to some finite time that has been set equal to zero without loss of 
generality; the corresponding configurations have well defined ultra-violet limit and they 
fully collapse to a point at t = 0. On the other hand, for c > 0, the scaling solutions are 
expanding as t runs from to oo and exhibit a well-defined infra-red limit. 

Irrespective of the uses and interpretation of the scaling configurations, their in- 
dependence follows by seeking solutions of equation (14.111) . Alternatively, the corre- 
sponding curves can be described using the parametric form of the mean curvature flow 

as 

Hh = cr (4.14) 

in terms of their position vector r = (x, y). As such, they may also be viewed as geodesies 
on the plane endowed with the metric 

ds 2 = e c{x2+y2 \dx 2 + dy 2 ) , (4.15) 

which is not Ricci flat but it is induced by the potential (14.121) according to the variational 
method of section 2.2 in the presence of dilaton. These are the scaling solutions that 
characterize the extrema of Huisken's entropy functional. In the physics literature, they 
first arose in the early work [H] and later in magneto-hydrodynamic models for the solar 
flares, [54] , where they were discussed in moderate detail. 

Another equivalent description is obtained by considering the mean curvature flow 
in polar coordinates, as in equation (13.16j) . with factorized t-dependence p(9(t),t) = 
R(9)^/2ct. Then, the slope of the curves depends only on 9, as (3(9), and satisfies the 
equation 

' -cR 2 {9) . (4.16) 



d9 

This formulation is advantageous for drawing the shape of the homothetic curves at fixed 
t. For this, let us assume without loss of generality that the curves are placed on the 
plane in a way so that (3(tt/2) = 0, meeting the y-axis perpendicularly at some point. 
Then, simple integration of equation (I4.16P yields 

(3(9) = c I"' 2 R 2 (9)d9 , (4.17) 
Je 

stating that the slope of any such curve is proportional to the area subtended by the 
corresponding radius vector as it moves away from its vertical reference position. Clearly, 
the curves are placed symmetrically about the y-axis, since (3 — ► —(3 when 9 — » 7r — 9. 
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Finally, note that the mean curvature of scaling solutions factorizes as H(/3, t) = 
H(P) I ^/2ct, where their dependence on the slope follows from equation (I3.13p . which 
now reads as 

HP + ™ + m-°- (4 - 18) 

This last equation turns out to be particularly advantageous for understanding the struc- 
ture of the homothetic solutions in detail as it provides an intuitive account for their 
classification depending on the sign of c. In the following, we study separately the self- 
shrinking and self-expanding solutions by stripping off their t-dependence and draw some 
characteristic figures that arise in each case. 

(i) Self-shrinkers (c < 0): It is apparent from equation (14. 18f) that solutions with 
constant mean curvature H(/3) can only exist for c < 0. They represent circles with 
radius R = 1/ \f—c for which H = sj—c. In this case, equation (I4.16P is trivially satisfied 
since (3 = 9 + tt/2 at all points of a circle. When time dependence is reinstated, the 
circles evolve by uniform contraction, as R(t) = y/2ct, until they collapse to a point; this 
special solution was studied in the context of boundary interactions in Ref . [26] , where it 
is referred to as circular brane model. Other solutions include rosette-like curves, which 
are symmetric about their maxima and minima, but they are not necessarily closed. 

There is a special class of solutions, however, which are closed rosettes with winding 
number p and q petals called Abresch-Langer curves r Pi9 , (58], [59]. Such curves are 
graphs of transcendental functions associated to any pair of relatively prime integers 
(p, q) so that 

The simplest one has characteristic integers (2, 3) and it is depicted in Fig. 5b next to 
the homothetically contracting round circle, whereas the next more complicated curve 
(3, 5) is depicted in Fig. 5c. Other examples correspond to the values (5, 8), (7, 10), (9, 14) 
(12, 17) and so on. 



(a) (b) 



Figure 5: Closed curves representing scaling solutions on the plane 
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The importance of these configurations stems from the fact that the mean curvature 
flow tends to evolve closed curves towards the scaling solutions with c < provided 
that the maximum curvature remains bounded as \H max (t)\ \r2ct < C < oo (type I 
singularities). Thus, embedded closed curves have the tendency to become circular as 
they shrink, whereas closed curves with self-intersections tend towards one of the special 
locally convex curves r Pj(? . The corresponding scaling solutions may be alternatively 
viewed as describing the asymptotic limit of the normalized mean curvature flow in the 
type I case. Furthermore, it can be shown under appropriate technical condition, that 
any closed n- dimensional hypersurfaces in R n+1 with non-negative mean curvature also 
evolves towards scaling solutions, [48J . Such solutions obey Hn = cf in all dimensions, 
after extracting their factorized t-dependence, as before, and they fall into three different 
classes: S n , S n ~ m x R m or T P:q x R n " 1 (see also Ref. [H]). Thus, the classification of 
scaling solutions on the plane for c < has more general value for the whole subject. 

We illuminate the presentation with a brief description of the transcendental nature 
of the closed Abresch-Langer curves r p (J , setting c = — 1 without loss of generality. First, 
note that equation ( 14.181) has a first integral 



E can be viewed as the energy of a point particle that moves with respect to an effective 
time (3 in a potential well V(H) having infinite height on both sides of the allowed range 
< H < oo. Since the minimum of the effective potential is reached at H = 1, in which 
case V(H) = 1/2, bounded motion with respect to (3 becomes possible for all E > 1/2, 
as in Fig. 6. 




(4.20) 



with integration constant E and 




(4.21) 



V(H) 



E 




2 



H_(E) 1 H+(E) H 

Figure 6: The effective potential for self-shrinkers 
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Periodic solutions H(/3) with fixed E have primitive period 



T(E) = 2 / . (4 22) 

J2(E-V(H)) 




where H±(E) are the roots of the equation = V^i?) that correspond to the two turning 
points of the bounded motion. The primitive period determines the minimum effective 
time that takes H to return back to its original value. This, however, does not necessarily 
mean that the curve itself has the same period in (3 for it may require several such turns, 
say q, for x(/3) and y{(3) to return back to themselves. Also, in other cases, the resulting 
curves may never close back to themselves, thus leading to rosette-type open shapes that 
wind indefinitely on the plane. 

Analysis of the problem shows that T(E) varies monotonically with E and it decreases 
from a/27t to n as the energy ranges inl/2<i?<oo. The upper bound of T(E) is easily 
established by considering H — 1+e with small but non- vanishing e. In this case equation 
(I4.20p takes the harmonic oscillator form 

(4.23) 

by expanding V(H) to quadratic order. Then, e{(3) is given by the trigonometric functions 
sjE - 1/2 sin(v / 2/?) or yj E - 1/2 cos(^2[3) with E - 1/2 small but strictly positive 
constant. These are periodic functions with period \/2tt that is insensitive to the value 
of E provided that E stays close to 1/2. The lower bound of T(E) can be established by 
asymptotic analysis that is rather involved and we refer the reader to the literature for 
the details, 158] , [59] . 



According to this result, closed curves on the plane with 1/2 < E < oo correspond 
to trajectories with primitive period 

T(E) = 2vr- , (4.24) 
Q 

where p and q are relatively prime integers subject to equation (14.19p above. This is a 
transcendental quantization condition for the parameter E showing that the emergence 
of closed rosettes on the plane is the exception rather than the rule. Other values of E 
result to rosette-type curves with infinite number of self-intersection points that never 
close back to themselves. Extending H{(3) periodically to [0, 2np\ yields the Abresch- 
Langer curves T p q with winding number p, as required, when the quantization condition 
( I4.24p is fulfilled. These curves have 2q critical points for their mean curvature reaches 
the minimal value H = 1 exactly twice within the primitive period T(E). They also 
appear to have q(p — 1) self-intersection points in general. For other values of E the 
homothetic solutions are open and can be formally thought to arise as limiting cases 
p — ► oo and q — > oo with infinite period. 

Finally, note that H(j3) follows from equation (14.201) by expressing (3 as indefinite 
integral of dH/ J2(E — V(H)), which obviously can not be written in terms of elementary 
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functions. The dependence of the position vector of such curves upon /?, as given by 
equation (13 . 1 21) . turns out to be transcendental. It is also clear in this context that an 
elementary solution arises when the effective point particle sits still at the minimum of 
the potential V(H) having E = 1/2. It corresponds to a round circle of unit radius so 
that its mean curvature is constant, H = 1, for all /?, and describes the only simple closed 
curve that shrinks by scaling in t. In that case the period of the curve is 2n and suffers 
discontinuous jump from the lower bound of T(E) when E > 1/2. 

(ii) Self-expanders (c > 0): This case allows only for hyperbola-like curves whose 
slope increases monotonically as one varies clockwise the polar coordinate 9. These 
curves are necessarily open with asymptotic lines placed symmetrically about the y- 
shown in Fig. 7. The need for asymptotic lines follows by inspection of the 
integral equation (14. 17ft . for, otherwise, the area subtended by the radius vector of the 
corresponding curve, as measured from the tip of the asymptotic wedge, will become 
unlimited in contradiction with the finite change of their slope. The curves can be 
thought as representing an intermediate stage for the decay of a wedge to straight line. 



Figure 7: Scaling solution on the plane representing the decay of a wedge 

The details of the solution can be investigated from the point of view of equation 
( I4.18p . setting c = 1 without loss of generality. As before, it has a first integral of the 
form (14.201) . with integration constant E, but with effective potential 



that differs from the c < case by a relative sign. As a result, V(H) is not bounded 
from below, for it is a monotonically increasing function ranging from — oo to oo when 
< H < oo. E can take any arbitrary real value in this case, and, clearly, V(H) can 
not support bounded motion with finite period; the corresponding scaling solution of the 
mean curvature flow is an open curve. For any given E, H(f3) follows, as before, by 



expressing (3 as indefinite integral of dH/^j2(E — V(H)), which in turn determines the 
form of the solution implicitly via transcendental functions. 

The mechanical analogue also helps to provide an intuitive explanation for the pres- 
ence of asymptotic lines and illustrate how scaling solutions can emerge from a wedge. 
H(f3) is the classical trajectory of a particle that rolls down the potential (14.251) having 
fixed energy E with respect to the effective time (3, as shown in Fig. 8. 





(4.25) 
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Figure 8: The effective potential for self-expanders 



The effective time that takes such a particle to go down the drain, when it starts with 
zero velocity, is 



where H (E) is the (single) root of the equation E = V(H) that specifies its position at 
some initial time /3 - Thus, as one transverses the corresponding curve on the plane, the 
slope changes by an overall finite amount A/3(E) that depends on E. The curve tends 
asymptotically to a straight line, with H = 0, which represents the universal attractor 
point of the potential. Actually, this effective particle motion traces only half of the 
corresponding curve on the plane, which can always be arranged to meet perpendicularly 
the y-axis by choosing /3 = 0; the mean curvature of the curve at the starting point is 
H (E), whereas the slope of the right asymptotic line is A(3(E) . The other half of the 
curve, together with the left asymptotic line, are placed symmetrically about the y-axis 
and can be obtained by simply reversing the direction of effective time. Then, tracing 
the complete curve from one asymptotic limit to the other amounts to shooting a particle 
up from the bottom of the potential and let it fall back to it after reaching a maximum 
height H (E) that depends on its energy. The two asymptotic lines meet at a point of the 
y-axis forming a wedge with opening angle n — 2A/3(E); this angle varies monotonically 
from Ti to as E ranges from — oo to oo. 

Particles with very low energy stay deep inside the throat of the potential, having 
H ~ everywhere, and they correspond to straight lines; in this case the wedge is wide 
open to 180° and the curve is lying horizontally on it. On the other hand, highly energetic 
particles have H (E) — > oo and the two sides of the wedge tend to collapse against each 
other. The same thing happens to the curve that folds up on the wedge, having infinite 
curvature at the vertex and practically zero everywhere else. For intermediate energies 
the typical scaling solution is a hyperbola-like curve with finite curvature everywhere. 
When the t-dependence is reinstated into the solutions, the coordinates x and y scale 
in the same way, as l/Vt, without affecting the angle of the asymptotic wedge y ~ \x\. 
The curves themselves appear as straight lines when one zooms closer and closer to 
them, whereas it becomes increasingly difficult to distinguish them from the surrounding 




(4.26) 
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asymptotic wedge when they are looked up from larger and larger distances away. 

Therefore, in this context, the scaling solution can be thought as the mean curvature 
analogue of the fundamental solution to the heat equation, whose initial configuration 
at t = is a delta function. It is known as Brakke's wedge in mathematics, |16j . |19j . 
but it also arose earlier in physics in Ref. [H] and [31]. The initial curvature singularity 
is fully dissipated after infinitely long time by flowing to the infra-red region. Clearly, 
it serves as model for studying tachyon condensation for intersecting branes and can be 
further used in connection with other worka^] on the subject (see, for instance, [50] and 
references therein). An analogous solution that describes the decay of a cone to the 
plane also exists for the Ricci flow, [61] (but see also [53]); it serves as model for studying 
tachyon condensation in closed string theory. Here, however, it is not possible to obtain 
the solution in closed form. 

4.5 Paper-clip model 

A genuine running solution of the mean curvature flow (13. 7p . with £ = 0, corresponds to 
the time dependent curve y = <p(x(t),t), 

e" 2t cosh.(vy) = cos(vx) . (4.27) 

with t running from — oo to some finite value T that has been chosen to be zero for 
convenience. The parameter v is free to take any arbitrary value. An equivalent form of 
the curve is 




(4.28) 



with two branches that are simply related to each other by y — > — y or equivalently by 
v — > —v. The variable x assumes values within the interval —7r/2v to ir/2v, but the 
precise range depends non-linearly on time t. 

The complete curve is closed since the two branches are glued symmetrically about 
the y = axis. For v = 0, the configuration reduces to a round circle that evolves by 
scaling of its radius, i.e., x 2 + y 2 = —2t, as can be seen by expanding cos(vx) and cosh(w?/) 
up to second order in their arguments. As such, it is common to the circular homothetic 
solution discussed earlier. For v ^ 0, however, the solution represents a convex curve in 
R 2 having oval (or paper-clip) shape at any given time t, as depicted in Fig.9; hence the 
name paper-clip model. 




y 



Figure 9: The paper-clip curve on the plane 
3 We thank Vassilis Niarchos for a discussion on this subject. 



32 



This solution arose independently in the physics and mathematics literature. It ap- 
peared as model for boundary interactions in Ref. [27], [2B], but it was also investigated 
much earlier in connection to magneto-hydrodynamics, [5l]. It has the special property 
that its extrinsic curvature H, when viewed as function of the slope (3 and the time t, 
satisfies the special ansatz 

H 2 {(3,t) = a(J3) + b{t) . (4.29) 
Substituting into equation (13.131) one finds the following system, 

(daV „ i a db , 2 , s 

y +4 °= 4c ' dt- 2b =- 2c > (4 - 3o) 

where c is a constant. The relevant solution in this class has 

2 2 

a (f3) = yCos(2/3) , b(t) = yCoth(-v 2 t) (4.31) 

for c = v 2 /2 and gives rise to the following dependence of the position vector upon (3, 

x{0) = - arcsin ^yfl — e 2v2t sin/3^ , y(/3) = — arcsinh C\J e~ 2vH — 1 cos/3^ , (4.32) 

up to translations, according to equation (I3.12p . Then, it is straightforward to verify 
that this coincides with the paper-clip curve ( 14.271) after eliminating the dependence on 
0. 

The corresponding arc-length, as measured from the tip of the paper-clip, follows by 
integrating dl = d(3 / H and equals to the incomplete elliptic integral of the first kind 

W = ~f 1 ^' i = ~m k) (4.33) 



in terms of the slope (3, with modulus k = v 1 — e 2v2t ; k varies from 1 to as t runs from 
— oo to 0. Thus, in general, we have the relation 

sin/3 = sn(vl/k; k) (4.34) 

in terms of the corresponding sine-amplitude Jacobi elliptic function. 

Clearly, when v — > one obtains the characteristic limit H 2 = —l/2t of a circular 
curve for all t. Also, when v 0, the ultra-violet limit t — > — oo of the paper-clip becomes 
asymptotic to the curve 

y = -log (2cos(ra)) — vt , (4.35) 

as viewed from the tip of the configuration associated to the y + branch. This is the 
hair-pin solution up to an irrelevant constant shift in t. In the ultra-violet region, the 
mean curvature of the paper-clip tends to the limit H 2 = v 2 cos 2 (3 with j3 ~ vx up to 
exponentially small corrections in time. A hair-pin facing in the opposite direction is also 
obtained in the ultra-violet region by viewing the curve from the other tip associated to 
the y- branch. 
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More generally, as t runs from — oo to 0, the paper-clip evolves by shrinking until it 
fully collapses to the point x — — yatt — and becomes extinct. Note that its size in 
the x-direction varies as 2v _1 arccos(exp(f 2 t)), and deminishes from n/v at t = — oo to 
zero at t = 0. Likewise, its size in the ^-direction varies as 2v _1 arccosh(exp(— v 2 t)), which 
also diminishes from infinite to zero length. Thus, as time goes on, the configuration 
becomes rounder and rounder by shrinking until it crunches to a point. Only when 
v = the two characteristic lengths of the configuration are equal and diminish evenly 
by preserving the circular shape of the corresponding solution. According to this, the 
paper-clip provides the mean curvature analogue of the sausage model encountered in 
the Ricci flow on5 2 , [H]. 

When the configuration is viewed from its "center of mass", and not from its tips, 
it looks as a "two-body" problem: two hair-pins with opposite orientation are glued 
together in their asymptotic region, y — 0, and move against each other until their tips 
merge. Their bound state is pictured schematically in Fig. 10 by putting together two 
periodic arrays of hair-pin and anti-hair-pin curves. The paper-clip is depicted by solid 
lines, whereas the dotted lines denote periodic repetition of the same configuration. 



i 
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Figure 10: Schematic gluing of hair-pins to form a paper-clip 

Note that a similar interpretation holds for the sausage model as a "two-body" prob- 
lem for Ricci solitons that merge together, [JT] , [53] . The metric of the sausage model can 
also be expressed in terms of the sine-amplitude Jacobi elliptic function, when written 
in proper coordinates, which is analogous to the elliptic dependence of the arc-length 
induced on the paper-clip curve. 

4.6 Oxlip model 

An open variant of the paper-clip model is obtained by considering the following curve 
on the plane, 

e v 'sinh(uy) = cos(vx) , (4.36) 
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for v 7^ 0. In this case, x ranges from —ir/v to ir/v, irrespective of t, and y can be either 
positive or negative. An equivalent description is given by the graph of the function 



y = -log ^cos(t>:r) + ycos 2 (vx) + e 2v2t j — vt . (4.37) 

For —ii/2v < x < ir/2v the variable y is positive, assuming that v > 0, whereas for 
—n/v < x < —n/2v and ir/2v < x < n/v the y coordinate is negative. The points 
(x = ±n/2v,y = 0) are inert under the flow, since these are the points of inflection 
where the extrinsic curvature vanishes. 

This particular solution appears to be new and it is quite interesting in many respects. 
When compared to the usual paper-clip model, the curve is formed by gluing together 
two periodic arrays of hair-pins facing in the opposite direction but with a relative shift 
in x equal to n/v , as in Fig. 11. Thus, it is a again a bound state problem of hair-pin and 
anti-hair-pin curves but of slightly different kind. The oxlip curve is depicted by the solid 
line, extending from —n/v to n/v, whereas the dotted lines denote periodic repetition of 
the same configuration. 




Figure 11: Schematic gluing of hair-pins to form an oxlip 

The remarkable feature of this solution is that the mean curvature flow exists for all 
time, as t ranges from — oo to oo. In the ultra-violet limit, t — > — oo, one readily gets 
an infinitely long hair-pin configuration as seen from its tip situated infinitely far away 
from the gluing region y = 0. On the other hand, in the infra-red limit t — > oo, the 
configuration tends towards the special curve y = with — n/v < x < n/v, which is a 
segment of a straight line. Thus, it appears that the hair-pin on the right is decaying 
while its two sides are getting squeezed against each other. Throughout the process, the 
two ends of the hair-pin stay firm at (x = ±n/2v,y = 0). The two half hair-pins on the 
left undergo similar decay until a line segment is finally formed in the infra-red limit. 
Three consecutive steps of the flow are depicted in Fig. 12. The same picture arises for 
v < 0, but with opposite orientation for the constituent hair-pins. 
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(a) (b) (c) 

Figure 12: Transition of a hair-pin to a line segment 



This solution satisfies the special ansatz H 2 ((3, t) = a(P) + b(t) for the mean curvature 
with 

2 2 

a(p) = ^-cos(2/5) , b(t) = ^-t&nh(-v 2 t) . (4.38) 

It is similar to the form of the paper-clip model, but b(t) is now given by the hyperbolic 
tangent rather than cotangent function; as a result, b(t) never blows up and the solution 
exists for all time. Then, x(/3) and y(P) can be calculated as for the paper-clip curve, and 
the same applies to the arc-length l(/3) which is also expressed in terms of the incomplete 
elliptic integral of the first kind. The only difference is the choice of modulus, which here 
is k = yl + e 2v2t , and varies from 1 to oo as t ranges from — oo to oo. It relates to the 
modulus k of the paper-clip curve by the transformation k = \Jl — k 2 in the common 
time interval (— oo, 0]. Standard identities among Jacobi elliptic functions show that the 
slope of the decaying hair-pin at different points of the curve is expressed in terms of its 
arc-length at any given instance of time by 

sin/5 = isn(vl; I) (4.39) 
k k 

with < \jk < 1. Clearly, the slope (3 vanishes everywhere in the infra-red limit where 
l/k = 0. 

Next, we will find, among other things, that the mean curvature flows for the paper- 
clip and oxlip models are naturally related to certain instability modes of the hair-pin 
configuration as one begins to move away from the ultra-violet region. 



5 Modes of instability and transitions on the plane 

In this section we perform a general analysis regarding the linearized (in) stability of 
special solutions that move by translations, rotations or scaling on the plane. In all 
cases, we work with the parabolic equation for the mean curvature and obtain an eigen- 
value problem of the form 

Ch{(3) = \h{(3) (5.1) 
governing small fluctuations about a given solution, 

H(0, t) = H (P) + h(P)esxp{Xt) . (5.2) 

Hq(P) denotes the mean curvature of the unperturbed curve in a reference frame that all 
time dependence has been accounted by the appropriate choice of vector field The 
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operator C is linear of second order, thus leading to an effective quantum mechanics 
problem in one dimension parametrized by the slope f3. The form of the corresponding 
potential depends on H {j3) and will be determined in all three cases. General conclusions 
about its spectrum will be drawn in each case separately. 

In general, eigen-states with positive A account for instabilities since the perturbations 
grow larger as time goes on. On the other hand, perturbations with negative A tend to 
dissipate exponentially fast and the configuration settles quickly back to its initial form, 
thus leading to stability. Usually, we refer to them as relevant and irrelevant deformations, 
respectively. Finally, if zero modes are present, they will only act on the parameter 
space of the given solution without affecting its time dependence. Our primary aim here 
is to identify potential instabilities about some given reference curves, which arise as 
fixed points of the mean curvature flow modulo translations, rotations or dilations, and 
associate transitions towards more stable configurations. 

5.1 Translating solution 

In this case, small fluctuations of the hair-pin take the form 

H(J3, t) = vcosp + h([3)e xt , (5.3) 

measuring its response against possible squeezing modes. Then, in the linearized approx- 
imation, equation (13.141) for the mean curvature becomes 

v 2 cos 2 (3 (fh + h)= Xh (5.4) 

since £ ■ h = vcos/3 and (d 2 /dj3 2 + 1)(£ • n) = 0. It turns out that the spectrum of A can 
be fully determined by simple transformation to an exactly solvable problem. 

For this, consider the change of variables 



z = log ^tan ^ + , ^(z) = h((3)Vcoshz (5.5) 

that transform the fluctuation equation (15. 4p into the following Schrodinger problem 

d 2 ^ 1 / 3 \ A _ , r n . 

+ 7 1 nr * = -V& ( 5 - 6 ) 



dz 2 ' 4 V cosh 2 z) v 2 

with — oo < z < oo as —tt/2 < (3 < tt/2. The variable z is obtained by integrating 
dz = d/3/ 'cos/3 = vd/3/H(/3) for the hair-pin curve and as such it coincides with its arc- 
length / (multiplied by v ) as measured from the tip; a useful relation here is cosh;? = 
1 / cos/?. Thus, we arrive at an eigen- value problem for a particle moving on the real line 
under the influence of a symmetric Rosen-Morse potential U(z) = W 2 (z) — W'(z) with 
superpotential 

W{z) = -tanhz (5.7) 
2 
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and energy E = —X/v 2 . 

This problem has been exactly solved in the literature, [62], in the context of super- 
symmetric quantum mechanics, and it was found that it admits only one normalizable 
eigen-state, 

*o(*) = ~7^= , (5-8) 
Vcoshz 

whose energy is zero. The ground state, which corresponds to h (P) = cos/5 under the 
change of variables (15.51) . does not induce a decay of the initial configuration but only 
amounts to shifting v by constant, as can be readily seen from equation (15.31) ; it is the 
expected behavior for a zero mode acting on the moduli v of the underlying hair-pin 
curve. There is also a continuum of scattering eigen-states with energies E > 1/4, which, 
however, have A < and correspond to stability modes of the problem as one flows away 
from the ultra-violet regime. According to this, the hair-pin (14. 3[) looks absolutely stable 
and acts as infra-red attractor for all hair-pin-like shapes that deviate from it by the 
appropriate perturbations. It is the expected behavior for a solitonic solution. 

There is an additional infinite set of discrete eigen-states in the Rosen-Morse poten- 
tial with E < 0, and hence A > 0, which may serve as instability modes of the hair-pin 
configuration. The existence of negative energy states in a problem of supersymmetric 
quantum mechanics looks strange at first sight, but it only implies that the corresponding 
Hamiltonian operator is not self-adjoint in the space of the corresponding wave functions, 
thus violating the lower energy bound. The reason is that such states are not normal- 
izable, since they blow up at the two ends of the hair-pin as z — > ±oo, and should be 
disregarded on normal grounds. This is indeed the case for the stability analysis of a 
single hair-pin, but as it turns out negative energy states play an interesting role in un- 
derstanding the linearized evolution of bound state problems of hair-pins, such as the 
paper-clip in section 4.5 and its open variant in section 4.6. It will be seen later that 
gluing two hair-pins with the opposite orientation amounts to canceling the divergent 
perturbations in their asymptotic region, thus giving rise to regular solutions. 

More precisely, it can be easily verified that that the Schrodinger problem ( 15.6P has 
eigen-states 

ty n ( z ) = i n +Wcoshz P*(ismhz) (5.9) 

with quantized energies E n = —n(n + 1) for all n — 1, 2, 3, • • -. Here, -P^(x) denote, up 
to normalization, the associated Legendre functions of x 

P^x) = VT^ ^[(1 - x 2 f] . (5.10) 

Despite appearances, all states ^ n (z) are real. Then, A assumes discrete values A n = 
v 2 n(n + 1), which are non-negative, and the corresponding modes are associated to in- 
stabilities of the hair-pin configuration. The ground state of the Rosen-Morse potential, 
^o(z) = l/Vcoshz, may also be appended to these formulae by extending their validity 
to n = 0. 
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The first non-normalizable state in this series has n — 1 with Ai = 2v 2 and ^i(z) = 
(coshz) 3 / 2 , which corresponds to hi{(3) = l/cos/3, up to normalization. A simple calcula- 
tion shows that this is precisely the perturbation driving the evolution of the paper-clip 
at the linearized level, as seen from one of its tips close to the ultra-violet region. Indeed, 
by expanding the mean curvature H(f3,t) of the paper-clip model one finds 



H(/3,t) = -^cos^) -coth(v 2 t) = vcos/3 + i^-f 2vH + 0{e AvH ) . (5.11) 



This is a valid expansion as long as t — > — oo and (3 parametrizes small deviations away 
from the tip of one of the two constituent hair-pins so that 1 / cos/3 remains bounded. It 
has precisely the form 



as noted above. Similar conclusions are drawn by expanding around the other tip of the 
curve. Likewise, for the oxlip model one finds 



which is valid in the same domain of parameters and corresponds to the same linearized 
perturbation (15. 12ft . 

This identification puts the negative energy states on firm basis and makes them 
physically relevant for the bound state problems under consideration. Of course, it is 
important to realize that the gluing conditions, which stick together the open ends of the 
constituent hair-pins, miraculously cancel the divergences of h\{0) on each component 
when (3 — > ±7r/2. The free ends of a hair-pin and an anti-hair-pin effectively attract each 
other, but this is a non-linear effect that can not be seen by expanding far away from 
their overlap. Only the exact solution, in either case, reveals the correct value of the 
extrinsic curvature at the connection points. Higher excited states n > 2 may also be 
used to describe potential decay channels for bound state problems of the hair-pin curves, 
since A n is strictly positive, but their form is more complicated; for instance, for n = 2, 
^2(2) = (coshz) 3 / 2 sinhz, which corresponds to h,2{f3) = sin/5/cos 2 /3 with A2 = 6t> 2 and 
so on. To the best of our knowledge, there is no exact description of the corresponding 
trajectories at the non-linear level, as in the case of the paper-clip and oxlip models. It 
is an interesting problem that deserves further investigation while searching for exact 
solutions of the mean curvature flow. 

Finally, we note for completeness that all positive energy states in the Rosen-Morse 
potential (I5.6P can be obtained from the discrete set of negative modes (15.91) by appro- 
priate continuation. More precisely, using the standard description of the associated 
Legendre functions in terms of hypergeometric functions, and replacing n by i k — 1/2 
in their arguments, so that E = —n(n + 1) becomes E = k 2 + 1/4, we obtain solutions 
with continuous spectrum E > 1/4 for all real values of k. If time were flowing in the 
opposite direction these would have been the instability modes of the problem. 



H(/3,t)~H {P) + h 1 (P)exp(\ 1 t) , 



(5.12) 





(5.13) 
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5.2 Scaling solutions 



Next, we investigate the existence of instability modes for perturbations of the scaling 
solutions. Naturally, there are two different cases here depending on the sign of the scaling 
parameter c. In either case it is appropriate to look at the linearized problem in terms of 
equation (13TT4"D for H((3,t) = H ((3) + h(f3)exp(Xt) with £-n = -c(x(J3) sin/? -y((3) cos/?) = 
-cS(0), in which case (d 2 /dp 2 + 1)(£ • h) = -c(S"((3) + S((3)) = -c/H((3). 

(i) Self-shrinkers (c < 0): For self-shrinking solutions the small fluctuation operator 
C reads 

£ c <o = tfo 2 (/3)J^ + #o(/?) + l (5-14) 

setting c = — 1 without loss of generality Focusing on closed curves, as we do in the 
sequel, amounts to solving the eigen-value problem Ch{f3) = \h(j3) on the space of 
periodic functions h(/3 + T) = h(f3) with period T = 2irp given in terms of the winding 
number p; note, however, that the periodicity on h{[3) is necessary but not sufficient 
condition to ensure closure of the resulting curve via equation (13.121) for x(j3) and y(/3), 
and it should be checked separately. Since the operator (15.141) is self-adjoint on the space 
of square-integrable functions L 2 (5' 1 ,(i/i) with measure dfi = dP/H^ft), it follows from 
the general theory, [63], that its spectrum is discrete 

Ao > Ai > A 2 > A 3 > A 4 > • ■ ■ (5.15) 

accumulating at — oo. The corresponding eigen-functions h n (f3) are orthonormal and 
although they can not be explicitly computed, due to the transcendental form of H Q (j3), 
they are bound to have a fixed number of nodes depending on n; in particular, h has 
no nodes, whereas hm-\ and h,2 n with n > 1 have exactly 2n zeros in S 1 within a period 
2-Kp. 

One easily sees in the present case that H (f3) is an eigen-f unction of the operator 
(I5.14p with eigen-value 2, which is necessarily the largest, i.e., A = 2, since h (j3) = 
Hq{(3) > 1 vanishes nowhere. Thus, there is at least one potential mode of instability 
modulo the question of keeping the deformed curve closed. There are additional modes 
of instability for the linearized perturbations of the Abresch-Langer curves T P)q with 
winding number p and q petals, which depend on q > 1. For this note that H' Q (/3) is also 
an eigen-function of the operator ( 15.141) but with zero eigen-value. Since it has 2q zeros, 
equal to the number of times Hq((3) reaches its minimum value 1 within a period 2irp, 
it follows from above that H' (f3) should be identified with either eigen-function h 2q -i 
or h 2q . Actually, it turns out that /?-2g-i(/3) = H' (f3) and zero is a simple eigen-value. 
The effect of the zero mode on the Abresch-Langer curves is to rotate them on the plane 
since Hq{(3 + e) = H Q (j3) + eH' (f3) to lowest order in e. Thus, in general, there are 
2q — l discrete modes of instability h , hi, • • • , h 2q -2 for all T p q curves; among these there 
are two eigen-functions, cos/? and sin/3, which both have A = 1. All other modes have 
negative eigen-values leading to exponentially damped perturbations in time. It turns 
out that the total number of instability modes that also preserve the closure of the curves 
is 2q — 3; further details can be found in Ref. [59] . 
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There is a more systematic way to examine the spectrum of small fluctuations using 



the arc-length I = z, for dz = df3/H(f3), and the variable ^(z) = h(j3) / y H Q (j3) , as for 
the translating solution. The effective quantum mechanics problem now reads 



The variable z ranges over a finite distance, equal to the total length of the T P:q curves, 
along which H remains positive ranging from H_(E) to H + (E). Thus, z can be regarded 
as periodic variable. Also, W(z) never becomes singular in this domain; it is equal to 
(one-half) the velocity of a particle moving in the potential well (14.211) . Thus, according 
to supersymmetric quantum mechanics, the energy spectrum is strictly positive, so that 
A < 2, supporting an infinite but discrete set of periodic solutions as discussed above. 
Of course, one may also have bands with continuum spectrum when more general Bloch- 
wave solutions are allowed to occur; these generalizations, however, do not yield periodic 
perturbations of the Abresch-Langer curves breaking their closure. 

The special case q = 1 is a round circle with H (/3) = 1 everywhere and its perturba- 
tions can be studied separately by solving the eigen- value problem (d 2 /d(3 2 + 2)h((3) = 
Xh(j3) on S 1 with period 2np. The periodic solutions are simply cos(n/3/p) and sin(n/3/p) 
with integer n. In either case, the eigen-values are A = 2 — (n/p) 2 , which are positive only 
for those n satisfying the inequality n 2 < 2p 2 . Clearly, there can be no zero mode for 
1 / \[2 is an irrational number; the absence of zero modes is also consistent with the fact 
that the position of the circle on the plane does not change by rotation, unlike the case 
of Abresch-Langer curves. For winding number p = 1, there are no unstable modes that 
preserve the closure of the curve. Indeed, although cos/? and sin/? have A = 1 and can be 
potential modes of instability, they are ruled out because they do not yield periodic x((3) 
and y(/3), respectively (see the remark towards the end of section 2.1). Thus, there is a 
unique simple homothetic closed curve on the plane which is absolutely stable against 
all perturbations; this is also consistent with the fact that a circle with winding number 
p = 1 always attracts locally convex simple closed curves under the flow, [46], [47]. For 
circles with higher winding number p > 1, it follows from above that the number of 
potentially unstable modes equals to the number of integers n, positive or negative, that 
satisfy the inequality \p/n\ > l/y/2. All these modes preserve the closure of the curve 
apart from two, cos/3 and sin/3, which have A = 1 for all p. 

Thus, as time goes on, one can envisage transitions from a circular configuration with 
winding number p towards an Abresch-Langer curve r Pi(? with 1/2 < p/q < 1/ \[2 followed 
by a transition towards a singular closed curved with q cusps, where the mean curvature 
blows up. The various stages of deformation are shown schematically in Fig. 13 for the 
simplest case (p, q) = (2,3), as in Ref. [58] . 





(5.16) 



where U (z) = W 2 



(z) + W'(z) with corresponding super-potential 




(5.17) 
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(a) (b) (c) 




(d) (e) (f) 

Figure 13: Transition from a double-folded circle to a singular curve 

In this case, turning on an instability mode, say h(/3) = sin(/3/2), for a circle with 
p — 2, which is pictured as two p — 1 copies sitting on top of each other, amounts to 
enlarging the mean curvature of the first copy (covered by < j3 < 2n) and decrease 
that of the second copy (covered 2ir < (3 < 4ir). It is easy to imagine that such uneven 
deformation, which keeps the curvature unchanged at the points (3 = 0, 2rr and An is 
pictured schematically by the closed curve Fig. 13b. Subsequent evolution towards the 
Abresch-Langer curve r2,3, as shown in Fig. 13c, is then due to non-linear effects. The 
instability modes of T2,3 are capable to deform it further towards the configurations 
shown in Fig. 13d and Fig.l3e until the singular curve shown in Fig.l3f is reached. None 
of the intermediate configurations correspond to scaling solutions apart from (a) and (c). 
Unfortunately, such transitions are not available in closed form; it is also an interesting 
problem for future work. 

The stability modes of the Abresch-Langer curves can be formally viewed as instability 
modes of the backward mean curvature flow. Thus, one may also envisage transitions 
from T p q to a circle with winding number p, as given schematically in the figure above 
by flowing continuously from Fig. 13c to Fig. 13b and finally to Fig. 13a. In this case, the 
periodic functions cos(q(3/p) and sm(q/3/p) with p/q < l/y/2 are stability modes of the 
circle with winding number p showing that the configuration Fig. 13a can act as attractor 
for the corresponding flow. It also explains from a more intuitive point of view the upper 
bound on p/q that defines the initial scaling solution F Ptq provide that such transitions 
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(among others) are indeed possible. In this respect, it can be shown that the mean 
curvature flow with initial data r = r p q + en, given in terms of the position vector of any 
Abresch-Langer curve T p q with |e| small, tends asymptotically to an p-fold circle when 
e > and to a singular closed curve with q cusps when e < 0, [54] . 

There are additional transitions one may envisage. One possibility is to have a decay 
channel for the singular curve Fig.l3f, which can be roughly thought as bound state of 
three wedges held together by appropriate gluing conditions, to a diminishing simple 
circle. Although there is no explicit solution of this kind, one may imagine that each 
corner will decay to a smooth curve coming out of the wedge and all three local solutions 
can be patched together to form a smooth closed curve that will eventually shrink to 
a point. Combining this transition with the one depicted in Fig. 13, we arrive at the 
reasonable conclusion that a configuration can change its winding number, e.g., from 
p — 2 to p — 1, when it passes through a singular shape. Another possibility based on 
the physical idea of tachyon condensation is that any curve with self-intersections, like 
the Abresch-Langer curves, will cut itself and follow the decay channel of two intersecting 
lines in the vicinity of each self-intersection pointQ. Recombination of the outgoing curves 
from each local wedge will result into a collection of circular branes that eventually 
shrink to points. This procedure resembles the construction of knot invariants out of 
planar closed curves, using appropriate cutting rules at the self-intersection points; see, 
for instance, [65] and references therein. These are depicted in Fig. 14 below. 




(a) (b) 
Figure 14: Cutting rules familiar from knot theory 

It may very well be that such formal connection holds the key for the systematic 
construction of the corresponding boundary states in quantum field theory. We plan to 
return to this problem elsewhere. 

(ii) Self-expanders (c > 0): Finally, we consider the case of self-expanding solutions 
for which the small fluctuation operator £ reads 

£ c>0 = + H 2 ({3) - 1 , (5.18) 

4 We thank Boris Pioline for a discussion on this point. 
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setting c = 1 for convenience. Note that the spectrum of £ c >o follows by subtracting two 
units from the spectrum of £ c <o, although the form of the eigen-functions h((3) is different 
for Hq{(3) is not the same. Also, in this case, the reference curves are open and one should 
impose a different set of boundary conditions on the corresponding eigen-functions. Here, 
Hq(P) is an eigen-function with A = —2, whereas H' Q ((3) is still an eigen-function with 
A = 0. This zero mode, as before, describes the freedom to rotate the curve on the 
plane and orient it differently with respect to the y-axis; the same rotation applies to its 
asymptotic lines. There are also two obvious modes with A = — 1 corresponding to the 
eigen-functions cos/3 and sin/3. 

In general we would like to solve the eigen- value problem (I5.18P with a prescribed set 
of boundary conditions. Here, (3 ranges in the interval [tt — A(3(E), A (3(E)] whose end 
points are given in terms of the slope (14.261) of the asymptotic lines to the scaling solution. 
If the perturbed curve is to remain asymptotic to the wedge, one has to consider solutions 
h{(3) that vanish at the end points of this interval. This is a well-defined bound state 
problem having discrete spectrum A < —2 that accumalates to — oo. The proof relies 
on the observation that Hq(/3), which is an eigen-state of the operator ( 15.181) with the 
above boundary conditions, vanishes nowhere else but at the end points of the interval; 
as such, it serves as the ground state of the problem. All higher excited states will have 
A < —2 and exhibit additional zeros at various intermediate points of the interval. The 
conclusion is that all perturbations of this kind correspond to stability modes. 

Instability modes may only arise if one alters the boundary conditions. Note that 
perturbations which can grow infinite large at the end points of the interval, while they 
remain bounded in all intermediate points, are capable to produce solutions with A > —2. 
For example, the zero mode H' Q {(3) is a monotonically decreasing function starting from 
oo at (3 = 7i — A(3{E) and ending to — oo at — A(3(E); it has only one zero at [3 — 
around which the perturbation stays small. It is easy to anticipate the existence of 
excited eigen-states with the same blow up behavior at the end points of the interval 
that exhibit more zeros at intermediate points. Although it is not possible to construct 
them explicitly, they are bound to have A > and lead to instabilities of the scaling 
solution. These positive modes are similar in nature to the instability modes of the 
hair-pin configuration that also blow up at the end points of the /3-interval. It will be 
interesting to examine them further and associate various decay channels towards some 
more stable configurations in "bound state" problems with curvature singularities, e.g., 
configurations of closed curves with cusps. 

A more systematic description is also provided here by supersymmetric quantum 
mechanics using the arc-length I = z that ranges from — oo to oo as one traces the curve 



from its far left to its far right asymptotic lines. Letting ^(z) = h(j3) / y H Q (f3) , as before, 
the corresponding linearized equation becomes 



where U(z) = W 2 (z) + W'(z) with super-potential given by the same functional form 





(5.19) 



44 



(I5.17p . This super-potential equals to (one-half) the velocity of a particle in the un- 
bounded potential well (I4.25P and can take any real value. Under usual boundary condi- 
tions at z — > ±00 the spectrum is non-negative leading to A < —2 by supersymmetry. 



5.3 Rotating solution 

Next, we consider small fluctuations around the static Yin- Yang curve with curvature 
H (/3), satisfying the defining relation (I4.10p . and substitute into equation (13.141) . Note 
that in this case £,-h = w(x(/3)cos/3 + y(/3)sin/3) = toS'{(3) and so simple calculation yields 
(d 2 /d(3 2 + h) = w(S"(J3) + S(J3))' = -uH'(/3)/H 2 ((3). Then, the corresponding 

linearized problem reads Ch(f3) = Xh(/3) with 

c = H ^ + 4 +Hm - 2 -W)- (5 - 20) 

Clearly, H' Q ((3) is a zero mode that corresponds to the freedom to perform rigid rotation 
of the Yin- Yang curve on the plane. 

In order to examine the spectrum of the operator L it is convenient, once again, to 
introduce the arc-length parameter I = z and change variable to 



This results into an effective Schrodinger problem for a particle in the state ^l{z) with 
energy E = —A moving in the following potential 

U(z) = W 2 (z) - W\z) - H 2 (z) , (5.22) 

where 

W(z) - i (J-*?®. + —) -- ( d 2M_ + (5 23) 

W[Z) ~2\H(z) dz + H(z))~2{ d{3 + H((3))- ^ 

The effective coordinate z ranges from —00 to 00, as one traces the complete Yin- Yang 
curve from one end to the other, with z = corresponding to the point of inflection 
located at its center. The curvature remains bounded everywhere and tends to zero far 
away from the center of the spiral. 

If the term H 2 (z) were not present, the potential U(z) would support a zero energy 
state as well as higher energy states by supersymmetric quantum mechanics. Note, 
however, that H 2 (z) is a positive definite term that vanishes asymptotically as z —>■ ±00 
and lowers the potential everywhere. Thus, the energy spectrum of U(z) is also lower 
having at least one negative energy state that arises by shifting the energy of the "would 
be supersymmetric" ground state. It also supports a zero energy state, as noted above, 
which arises by shifting the energy of some otherwise excited supersymmetric state; 
furthermore, it supports other positive energy states. Although it is difficult to determine 
the exact form of the shifted energy eigen-states, it is clear that the presence of negative 
modes, which have A > 0, will lead to instabilities of the Yin- Yang curve. 
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6 Mean curvature flow on two-dimensional surfaces 



In this section we extend the previous discussion to two-dimensional curved spaces, so 
that the mean curvature flow is naturally combined with the Ricci flow. We will present 
examples of curves embedded in conformal backgrounds, such as the Ricci soliton, as well 
as examples of deforming curves embedded in spaces with deforming metrics. Among 
them there are mini-superspace models that capture the competition of shrinking curves 
on shrinking backgrounds. Closed curves may fully collapse before the metric reaches a 
singularity or they may collapse simultaneously to a point. 



6.1 General aspects 

Let us consider a general two-dimensional ambient space M. whose metric is expressed 
in conformally flat form 

ds 2 M = ^(dx 2 + dy 2 ). (6.1) 

The conformal factor Q depends on the coordinates x and y and it may also depend on 
the renormalization group time t when the background deforms under the Ricci flow. In 
this space we consider embedded curves (x(s),y(s)) whose induced line element is 

-•-s ((£)'+(£))*■■ (6 - 2) 

According to the definitions given in appendix A, the unit normal vector to each point 
of these curves is 

n = . = (-<p'(x), 1) , (6.3) 

Ji + ^Hx) 

and the mean curvature is 

H = , + 1 (d y Q - tf/(x)d x Q) , (6.4) 

(Jl + ip' 2 {x)\ 2Vn^i + v > 2 (x) 

which generalize the corresponding expressions for planar curves used in section 3. The 
notation y = (p(x) is used here, as before, to express the embedded curves in graph form. 

The mean curvature flow in M. is formulated, as usual, by computing the deformations 
of the coordinate functions x(s,t) and y(s,t) driven by the mean curvature vector. The 
result is better described by the deformation of graphs y = ip(x(t),t), which turn out to 
satisfy equation 

% = i+yi(!) + \ ^ - ^ x ) d * Q ) - kv + ^'^ kx ■ ( 6 - s ) 

Here, we have properly included the effect of reparametrizations generated by a Killing 
vector field k (if it is at all present) along the flow. This equation should be combined 
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with the Ricci flow of the target space metric, which reads in terms of f2 

dQ- 1 1 



Of 



o(3 + ^)iog(n -1 ) (6-6) 



when there is no dilaton present in the model. Thus, the combined system of Ricci 
and mean curvature flows for embedded curves in arbitrary two-dimensional Riemannian 
spaces are described by equations (16.5p and ( 16. 6ft for the two unknown functions Q and 
ip, which in general depend on t. 

In the presence of a dilaton field $, the conformal factor should satisfy the following 
constraints 

d s (nd y $) + d y (nd x $) = o , d x (nd s $) = d y (nd y $) , (6.7) 

for, otherwise, the different components of the Ricci flow equations for the metric are 
not compatible. In view of the applications that will be considered later, let us assume 
that both $ and Q are independent of the coordinate x so that the model exhibits an 
isometry with Killing vector field d/dx. Then, the constraints above are automatically 
satisfied provided that Qd y Q is independent of y. This term is actually constant, and not 
function of t, for otherwise the ^/-derivative of the dilaton flow will be incompatible with 
the Ricci flow. This constant is of order 1/a' and may be chosen so that 

fld v $ = — (6.8) 
a' 

when the dilaton field is non-trivial; otherwise it is zero. The normalization ensures that 
fixed points of the Ricci flow are also fixed points of the dilaton flow accounting for the 
balance between the field dependent and central terms of /?($) in two dimensions. Thus, 
the two fields <3> and Q are taken to satisfy equation (16.81) for all t, and, from now on, we 
set for convenience a' = 4. 

According to this, the Ricci flow for Q(y,t) is given by the evolution 

^n-^^togCn- 1 )-^^- 1 ), (6.9) 

whereas the dilaton $(y,t) follows by integration of equation (16.81) . Also, the equation 
for the mean curvature flow, expressed in terms of the graph (p(x(t),t) for embedded 
curves, takes the following form 

accounting also for the dilaton and the effect of reparametrizations generated by the 
Killing vector field (k x , k y ) along the flow. 

It is convenient sometimes, when the ambient space is surface of revolution, to use 
proper coordinates (r, 9) so that the metric takes the form 



ds 2 M = A 2 (t) (dr 2 + f{r, t)d9 2 ) (6.11) 
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with scale factor A(t) and profile function f(r,t). These coordinates are best suited for 
drawing the surface and the curves embedded in it. The corresponding Killing vector 
field is (k r ,k e ) = (0, o>) and it is associated to arbitrary angular velocity u. Assuming 
that the dilaton is independent of 9, the Ricci flow equations become 

f)A 2 f" ftf 1 

dt f dt A 2KJ J { ' 

where prime denotes the derivative with respect to the radial coordinate. The dilaton 
flow accompanying them is 

setting, once again, a' = 4. Also, the mean curvature flow for deforming curves can be 
formulated in this frame by considering r = p(9(t),t). Explicit calculation shows that it 
takes the following form 

^ = 1 (*p _ m ( f y 0) + 2 (?py\\ + + 

dt A 2 {p{p) + {d P /de) 2 ) yoe 2 f{ P ) y J Kl ' \de)j) a 2 yp> ae ' 

(6.14) 

adding also the contribution of the dilaton and the effect of possible reparametrizations 
generated by uniform rotation along the flow. 

Rotating solitons exist on all surfaces of revolution, thus generalizing the planar Yin- 
Yang curve. They correspond to fixed points of equation f)6.14p . with w ^ 0, satisfying 
the ordinary non-linear differential equation for p(6), 

^ arctan [mm) + !{pmp) " / ' ( " ) + ^ 2/(p) ^ = • (6 ' 15) 

This equation can not be easily solved in closed form. Even in the simplest case of 
planar rotating solitons, for which equation (I6.15P can be integrated once, the solution 
is only given implicitly. On curved surfaces, the background fields undergo continuous 
deformations by Ricci flow, and, therefore, rotating solitons correspond to curves solving 
equation (16. 15ft at each given instance of time. Some non-trivial examples will be studied 
later. 

Another class of special solutions arise on all surfaces of revolution when p is indepen- 
dent of the angular variable 9. They correspond to circular curves that can roll on the 
surface while remaining symmetric about the principal axis. In this case, p only depends 
on t and satisfies the equation 



dP 1 ^ fW 



dt A 2 V v " f(p) 



(6.16) 



The solutions generalize the uniformly contracting circular planar curves to curved back- 
grounds. However, their time evolution depends on the geometry of the surface. We will 
also see examples of this later. 
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6.2 Branes on constant curvature surfaces 

Trivial dilaton in proper coordinates implies that the profile function f(r) is independent 
of t, so that all time dependence is fully encoded into the scale factor A(t). In this case, 
f"(r)/f(r) is constant that can be normalized to —1, +1 or without loss of generality. 
These are precisely constant curvature metrics on Ai with positive, negative or zero 
curvature, respectively, for which A 2 (t) turns out to be — t, t or 1, up to normalization. 
Here, we will examine some simple solutions of the mean curvature flow on the sphere 
with f(r) = sinr and on the one-sheeted hyperboloid with f(r) = coshr. In the former 
case, r is an angular variable ranging from to it as one moves from one pole of the 
sphere to the other, whereas in the latter r ranges over the entire real line. 

The simplest running solutions correspond to circular curves satisfying equation (16.161) . 
On the sphere one gets the solution 

|cosp(f)| = -j , (6.17) 

where to is a non-negative integration constant. Thus, as time flows from — oo to the 
background is a uniformly contracting sphere. The circular brane appears to come from 
the equator and slip off the side, on either hemisphere, until it fully collapses to a point 
at t — —to before the big crunch. There is also the special solution p = it/ 2, which 
corresponds to t = and represents a great circle that follows the collapse of the sphere 
all the way to a point. The corresponding solution on the hyperboloid reads 

|sinhp(t)| = |, (6.18) 

where the integration constant t is again non-negative. In this time flows from 

to oo, the background evolves uniformly by lowering its curvature. The circular brane 
appears to come from the asymptotic region of the hyperboloid, which corresponds to 
p — > ±oo, and stabilizes to a circle at p = in the infra-red limit. As before, there is 
also the special solution p = for to — that does not roll at all on either side. 

Another interesting problem is the construction of rotating solutions on constant 
curvature spaces. The problem has already been investigated, to some extend, in Ref. 
[66] under the unnatural condition that the metric does not Ricci flow. These results, 
however, can be easily generalized to uniformly varying backgrounds without much effort. 
One way is to rescale the metric and redefine time so that constant curvature metrics 
appear as fixed points of the normalized Ricci flow. The mean curvature flow should 
be modified accordingly when expressed in the new variables. Another way is to use 
equation (16 . 15[) for rotating solitons, as it stands, and define an effective angular velocity, 

u cS = uA 2 . (6.19) 

Then, the shape of the resulting curves is identical to those drawn in Ref. [66] at any 
given instance of time. As time flows, u> e g changes; on the sphere it diminishes from oo 
to 0, whereas on the hyperboloid it increases from to oo. Following the analysis of Ref. 
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[66] , we give a schematic representation of the rotating solitons at some intermediate 
time. 

On the sphere, the curve oscillates around the equator and keeps coming closer to it 
as it winds. In the process it also keeps crossing itself. This behavior is best seen by the 
numerical plot r(6) shown in Fig. 15 below. 








Figure 15: The graph of a rotating soliton on the sphere 

On the hyperboloid, the curve resembles the shape of the planar Yin- Yang curve, 
which is now stretched on from — oo to oo along the symmetry axis. The corresponding 
solution is depicted in Fig. 16 below. 




Figure 16: A rotating soliton on the one-sheeted hyperboloid 



6.3 Branes on Ricci solitons 

The simplest static example of curved ambient space is provided by the following choice 
of metric and dilaton fields, in conformally flat frame, 

n(y) = 1 + e 2y , <%) = -Uog (l + e-^) . (6.20) 
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The coordinate y takes all values from — oo to oo, whereas x is an angular variable ranging 
from to 2ir; in this frame, the space is conformally equivalent to the cylinder^. The 
configuration corresponds to the well known Ricci soliton associated to non-trivial fixed 
point of the Ricci flow in the presence of dilaton, [5T] (but see also [8]). In the physics 
literature it serves as model for string propagation on a two-dimensional Euclidean black 
hole background, [55]. The mean curvature flow will be studied in this space using 
equation (I6.10P in the system of conformally flat coordinates. Since d y Q = 2Q — 2, the 
mean curvature flow simplifies to 

+ fcVfa) - k ' y • ( 6 -21) 

This background exhibits a rotational isometry generated by d/dx, so that k x = u and 
k y = with angular velocity u of either sign. Clearly, there is no other isometry. 

The fixed points of the mean curvature flow correspond to curves described by the 
equation 

ip"(x) 

1 I ,2 ; +1 = 0, 6.22 

1 + (f {x) 

without taking into account reparametrizations generated by the Killing vector field k. 
It is identical to the hair-pin equation on the plane, with parameter v — — 1, and the 
general solution is 

e y-y = cos ( x _ Xq ) ; ( 6>2 3) 

allowing also for the possibility to shift the coordinates by constant (x ,y ). Thus, 
we recover the standard Dl-brane on the Ricci soliton that exists by itself, without 
need for reparametrizations, thanks to the special form of the supporting dilaton field. 
This coincidence makes the hair-pin a rather special configuration with infinitely many 
symmetries inherited by the Euclidean black hole geometry. 

The same solution can be alternatively described using proper coordinates (r, 8) in 
target space. In this frame, the Ricci soliton corresponds to the choice of profile and 
dilaton functions 

f(r) = tanhr , $(r) = — log(coshr) , (6.24) 

whereas the overall scale factor A is a constant set equal to 1. Then, the hair-pin on the 
Ricci soliton is the curve r = p(0), 

sinhpo 



sinhp 



cos(# - Oo) , (6.25) 



which provides the static solution of equation (16.141) with uj = 0. Of course, the two 
descriptions are related to each other by the coordinate transformation 

sinhp = e~ y , 6 = x . (6.26) 



An alternative description exists by changing frame to X ± iY — exp(— y ± ix) that maps the 
cylinder to the plane. Then, the metric is conformally equivalent to the plane with metric ds 2 = 
(dX 2 + dY 2 )/{\ + X 2 + Y 2 ), whereas the dilaton is -2$ = log(l + X 2 + Y 2 ). 
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The tip of the cigar corresponds to y — oo and its asymptotic region to y = — oo. 

The parameters po and 8q are integration constants, but, clearly, p > po for the 
solution to make sense. Thus, po determines the position of the tip of the hair-pin 
relative to the tip of the cigar, which is located at p = 0. By the same token, y < y . 
The other parameter, 8 Q = xq, measures the rotation angle of the cigar relative to a given 
position about its axis and can be set equal to zero for all practical purposes. In the 
asymptotic region, p — > oo, the cigar looks like a cylinder and the hair-pin reduces to 
a pair of diametrically opposite parallel lines placed on it, since 6 — 9 — > ±7r/2. As p 
decreases, these two lines bend towards each other and meet smoothly at p = po- When 
Po = 0, the curve passes through the origin and its two legs are diametrically opposite for 
all p. The solution represents an open curve sitting still on the surface of a semi-infinite 
long cigar, as shown in Fig. 17 below. 



We will see later that this particular configuration can be used as component for 
constructing the analogue of a paper-clip on axially symmetric evolving backgrounds 
with spherical topology. 

Rotating solutions on the two-dimensional black hole geometry correspond to gener- 
alized fixed points of the mean curvature flow with uj ^ 0. Explicit calculation shows 
that the resulting equation can be integrated once, as on the plane, leading to 



The same analysis can be performed in proper coordinates, where the rotating soliton 
satisfies the equation 



To compare this curve to the planar Yin- Yang spiral, it is necessary to use a common 
frame. Thinking of the plane as being conformally equivalent to a cylinder (x, y) with 
periodic variable x and metric ds 2 = (dx 2 + dy 2 )/e 2y , the planar Yin- Yang curve satisfies 
the equation 




Figure 17: A hair-pin curve supported by the cigar background 




(6.27) 




(6.28) 




(6.29) 
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The two equations match on the side y — > oo, in which case the rotating soliton becomes 
independent of uj and approximates the hair-pin curve arctan<//(:r) + x = 0, i.e., exp(y — 
y ) = cosx. Note that the tip of the hair-pin should be placed very far away, i.e., 
y — > oo, for otherwise the approximation would have not been valid. According to this, 
the rotating solution on the cigar starts from its tip, as expected by symmetry. 

The identification of the two curves is also valid relatively close to infinity, since 
log(l + e~ 2y ) ~ e~ 2y , but it breaks down to order 0(e~ 4y ). Thus, close to the tip of the 
cigar, the equation becomes approximately 

u e~ 2y 

(f'(x) = -tana; + — — — , (6.30) 
2 cos z x 

dropping all terms of order C(e~ 4?; ), and it is solved by y = <p(x), 

e 2y = e 2?/o cos 2 x + ^( tanx + sin2x ) . (6.31) 

Setting 5y = y — yo and expanding the trigonometric functions around the tip, where 
x = 0, we obtain to first order 

Sy = ^ e - 2yo 5x . (6.32) 

This shows the tendency of the curve to twist as it moves away from the tip. 

On the other hand, close to the asymptotic region of the cigar, y — > — oo, the rotating 
soliton is described approximately by equation 

<p'(x) = taa.(uiy — x) . (6.33) 

It is solved exactly by the following expression 

y + ujx = log ^sm(uy — x) cos(u;?/ — x)j (6.34) 

up to an irrelevant integration constant. Then, it becomes clear that the asymptotic 
dependence of the curve y = f(x) is 

ujy — x = arctan— (mod 7r) , as y — > — oo . (6.35) 

UJ 

For, otherwise, the right-hand side of equation (I6.34p can not match the infinity appearing 
on the left-hand side. This assertion can also be verified by direct substitution of (16.351) 
into equation (I6.33p . When uj = 0, the asymptotic relation (16.35)) reproduces the well 
known asymptotic description of the hair-pin as two parallel lines with x = ir /2 (mod7r). 
When w ^ 0, the hair-pin ends up in a double helix whose components are tt/uj apart 
from each other. The sign of uj determines the handedness of the helix. 

The corresponding curve is centered at the tip of the cigar and winds around the 
asymptotic cylinder as shown in Fig. 18. The structure of the curve is more complicated 
in the middle region and thorough numerical analysis is required to draw its shape. 
Certainly, it is quite different from the planar Yin- Yang curve. 
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Figure 18: The end-point structure of a rotating soliton on the cigar 

Finally, we consider running solutions on the cigar that represent circular branes 
placed perpendicularly to the axis of symmetry. In this case p depends solely on t and 
simple integration of equation (16 . 16[) yields 

coshp = e' "' (6.36) 

when A — 1. These curves originate from the asymptotic region of the cigar at t = — oo 
and move towards the tip until they fully collapse to the point p = at some finite time 
t = to- They can be thought as analogue of the uniformly shrinking circles on the plane, 
although the radial dependence on t is different now. It is natural to expect that all 
closed curves winding around the cigar will contract to a point at its tip irrespective of 
initial conditions. It should be analogous to the well known fact on the plane that all 
closed curves shrink to a point at finite time. 



6.4 Branes on a sausage 

The sausage model was introduced in the physics literature more than a decade ago, [IT] , 
and describes axially symmetric deformations of the sphere by Ricci flow. It is defined 
by the following ansatz 

n{y, t) = a(t) + 6(t)cosh2y (6.37) 

in a system of conformally flat coordinates (x,y), whereas $(?/) = 0. The coordinate 
y can take all values on the real line and x is taken to be periodic ranging from to 
2tt. This ansatz yields a consistent truncation of the Ricci flow to a simpler system of 
ordinary differential equations for the two moduli a(t) and b(t), 

da , db , 

— = 2b 2 , - = 2ab , (6.38) 

which can be easily solved as 

a(t) = 7 coth(-2 7 t) , b(t) = . , 7 _ • (6.39) 

smh(— Z'yt) 

The integration constant 7 is assumed to be non-negative and determines the ultra-violet 
limit of the configuration. Indeed, as t — > —00, Q tends to a constant value, 7, and the 
space looks like an infinitely long cylinder of radius 1 / y/j. Then, as t increases, the 
configuration looks like ge that evolves by becoming shorter and rounder until 

it fully collapses to a point at some finite time t — 0. When 7 = 0, the trajectory 
corresponds to a uniformly contracting round sphere of radius 
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On this two-parameter space we will examine the mean curvature flow for embedded 
curves, y = (p(x(t),t), 

^ = — — r [ (cosh(-2 7 t) + cosh2y) - - + sinh2y ) (6.40) 

dt sinh(-2 7 t) V v v ' J y 'l + V ' 2 (x) y j V ; 

and construct special solutions that generalize the paper-clip and the oxlip planar curves 
to the sausage. They provide explicit realizations of the curve shortening problem on 
a deforming background and, in this respect, there are two natural scale parameters at 
work. There is the time at which the background is fully collapsed to a point, taken here 
to occur at t — 0, and the time at which the closed curves collapse to a point, denoted 
by t m the sequel. Clearly, t < and in most cases one expects, based on intuition, 
that the curve will become singular before the background. The paper-clip and other 
examples that will be considered later have t < 0, but for the oxlip model we find t = 0. 
The solutions simplify considerably when 7 = 0, in which case they describe deforming 
curves on a uniformly contracting sphere. 

It can be easily verified that equation (I6.40p admits the following simple solution 

KCOsh V ( \ (p. a~\\ 

cos(a; — x ) , (6.41) 



sinh(— jt) 

where k is a positive integration constant and xq represents the freedom to rotate the 
sausage by an arbitrary angle about its axis of symmetry. This solution describes a closed 
curve on the sausage that can be thought as superposition of two hair-pin solutions. As 
such, it provides the analogue of the paper-clip curve on the sausage. One way to see 
this is by rewriting the sausage model metric as 



2 1 / I I 

ds 1 



M 



ry \ I _|_ e 2y+2 7 t I e 2y-2",t 



) (dx 2 + dy 2 ) . (6.42) 



It represents the bound state of two cigars glued together in their asymptotic region that 
begin eating each other as time flows. The corresponding dilaton fields also come with 
opposite signs and cancel each other, up to spatially independent terms. Then, the closed 
curve (I6.4ip can be interpreted as two hair-pins placed appropriately against each other 
on the constituents cigars. The decomposition of the curve on very long sausages reads, 
in particular, 

[e y+ ^ + e- y+7t ) = cos(x - x ) + C(e 27 *) . (6.43) 



The gluing of the individual components is performed in a certain way, as depicted 
schematically in Fig. 19 below. Another type of gluing condition will be considered later 
that is reminiscent of the oxlip planar curve. 



Figure 19: The paper-clip as composition of two hair-pin curves on cigars 
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The picture becomes more clear in proper coordinates, which are best suited for 
drawing pictures. Introducing the change of variables 



sn(r; k) = tanhy , 9 = x , k — tanh(— jt) , (6.44) 
the sausage model metric assumes the following form 

ds 2 M = - (dr 2 + sn 2 (r; k)d6 2 ) (6.45) 

written in terms of the Jacobi elliptic function sn(r; k) with modulus k. One tip of the 
sausage is located at r = and the other at r = 2K(k) given in terms of the complete 
elliptic integral of the first kind. There is also a dilaton field 

$'(r) = rk' 2 — E(r, k) (6.46) 

whose form is determined by consistency of the Ricci flow in proper coordinates. Here, 
E(r, k) denotes the incomplete elliptic integral of the second kind; $(r, t) can be expressed 
in terms of Jacobi's theta function by simple integration. 

The ultra-violet limit of the ambient space corresponds to A; = 1, in which case 
the length of the sausage becomes infinite, whereas k = corresponds to t = and 
the configuration fully collapses to a point. In the ultra-violet limit, the sausage looks 
like a cylinder from its middle point y = 0. In proper coordinates, however, one sees 
an infinitely long sausage from one of its tips and the structure looks identical to the 
two-dimensional cigar, 

ds 2 M ~ -(dr 2 + tanhVdfl 2 ) , (6.47) 

7 

since sn(r; 1) = tanhr. The radius of the circle in the asymptotic region r — > oo is 1/^/7, 
as required. The picture is alike from the other tip, since sn 2 (r + 2K(k); k) = sn 2 (r; k). 
Thus, the interpretation of the sausage as bound state of two Euclidean black holes 
becomes rather precise. 

The paper-clip on the sausage takes the following form in proper coordinates, 

Kdn(r;k) 

k^k) = C ° S{6 " ° o) (6 ' 48) 
using Jacobi elliptic functions. In the ultra-violet limit k = 1 one has 

dn(r; 1) 1 



sn(r: 1) sinhr 



(6.49) 



and so one recovers the hair-pin curve (16.251) on one of the two constituent cigars with 
parameter k = sinhro- From the other tip of the sausage one sees a second hair-pin 
placed symmetrically with respect to its center, as in a mirror; on the second hair-pin 
one makes the identification k = — sinhr because the function dn(r; fc)/sn(r; k) flips sign 
when its argument is shifted by 2K(k). Then, the two pieces are glued together in the 
central region. The open ends of each hair-pin, which are diametrically opposite in the 
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ultra-violet limit, are joined smoothly to form an infinitely long paper-clip on an infinitely 
long sausage. As time goes on, the curve tends to slip off the side until it collapses to a 
point before the sausage shrinks to zero size. The collapse of the curve occurs at time to 



sinh(— t^o ) = k , i.e., k = , = . (6.50) 

sfl + 1/ 'k 2 

The parameters 7 and k are, in general, independent. Fig. 20 below depicts the 
evolution pattern of the paper-clip on a sausage with 7 7^ 0. The limiting case of a 
paper-clip on the uniformly contracting sphere appears to be singular because 7 = 0. It 
can only be accommodated by considering a correlated limit of parameters so that /c/7 
remains fixed to a constant. 





Figure 20: Schematic evolution of a paper-clip on the sausage 



An interesting variant of the solution above is provided by the following curve on the 
sausage, 

sinh?/ 

cos(x — x ) . (6.51) 



cosh(— 7^) 

As can be easily checked, it satisfies the mean curvature flow (16.401) . When expressed in 
proper coordinates it takes the form 



cn(r; k) 
sn(r; k) 



cos 



(6.52) 



In the ultra-violet limit, the curve looks similar to the one considered before, because 
cn(r; 1) = dn(r; 1). However, the picture is slightly different from the other tip of the 
sausage because both functions sn(r; k) and cn(r; k) flip sign when their argument is 
shifted by 2K(k). As a result, the relative orientation of the two hair-pins changes and 
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one gets k = sinhp on both sides. Thus, the solution represents a closed curve on the 
sausage, which is formed by putting together two hair-pins on the constituent cigars, as 
before, but this time the gluing prescription is different. In particular, one of the two 
cigars should be rotated by an angle 7r before gluing it to the other. This operation does 
not alter the metric of the ambient space, but affects the curves embedded in it. This is 
also clearly seen by comparing equation ( 16.5ip on very long sausages, 



to the analogous expression (16.431) for the paper-clip curve. The flip of the relative sign is 
attributed to the rotation of the second component. The resulting configuration, before 
and after the twist, is depicted in Fig.21 below. 



This solution provides the analogue of the oxlip curve on the sausage. The analogy is 
revealed by comparing equation (16.511) in conformally flat frame to the oxlip solution on 
the plane. They both share the characteristic sinhy dependence as opposed to the coshy 
dependence of the paper-clip curves. The oxlip solution exists for all time as long as the 
ambient space is regular. Thus, on the plane, it exists for infinitely long time, whereas on 
the sausage it shrinks together with the space until they become singular simultaneously 
at t — 0. This behavior is easily understood on intuitive grounds: the curve goes around 
the two tips of the sausage on opposite sides and cannot slip off to one side. As time goes 
on, the curve has the tendency to become shorter by moving evenly towards the center 
of the sausage, which also gets shorter and rounder, until the big crunch. The topology 
of the ambient space affects only the shape of the oxlip curve, which is closed on the 
sausage and open on the plane. 

7 Mean curvature flow in three dimensions 

In this section we consider the mean curvature flow of two-dimensional surfaces embedded 
in three dimensions. We will only examine the case of flat ambient space, R 3 , which 
already poses a non-trivial problem. Branes in curved ambient spaces are much more 
difficult to study since the mean curvature flow should be combined with the Ricci flow, 
which is very complex problem. We will avoid such unnecessary complications and also 
ignore the effect of fluxes that can be turned on in three (and higher) dimensional spaces. 




(6.53) 




Figure 21: Cross-joining two hair-pin curves on cigars 
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The role of fluxes in the boundary renormalization group equations of sigma models is 
fully captured by the Dirac-Born-Infeld that will be further studied elsewhere. 



7.1 General aspects of evolving branes in R 3 

The mean curvature flow of surfaces in R 3 is obtained by considering the general embed- 
ding equation 

X = X(s,u;t), Y = Y(s,u;t), Z = Z(s,u;t) (7.1) 

associated to two parameters s and u and the renormalization group time t. It is 
convenient, where appropriate, to think of the surface as graph of a function Z = 
(p(X(t),Y(t);t) that evolves in time. 

Using the formulae given in appendix A, it turns out that the mean curvature of the 
surface is 

H = (1 + (dy V f) d\ V + (1 + {d xV f) %r<p - 2{d X <p){dY<p){d X dY<p) 

( V /l + (^) 2 + (^) 2 ) 3 

whereas the inward unit normal vector is 

n = 1 {-dx¥, - dy(p, 1) • (7.3) 

'l + (d xV ) 2 + (dy^Y 



Then, the mean curvature flow takes the form 

(1 + (dyipf) & x ip + (1 + [dx^f) &YV ~ 2(dxV)(d YV )(dxd Y cp) 

dt " 1 + (dxv) 2 + (d Y cpy 

H X dxV + edytp - e (7.4) 

by also adding the effect of arbitrary reparametrizations generated by a vector field 
(£ x ,£ y ,£ z ), if appropriate. 

The induced metric on the two-dimensional surface evolves according to the equations 
^g AB = -2HK AB , ^g AB = 2HK AB , (7.5) 

so that 

-v/det^ = -H 2 Jfetg . (7.6) 
Also, the second fundamental form of the surface follows the evolution 

jK AB = g CD V c V D K AB - 2H(K 2 ) AB + (TrK 2 )K AB , (7.7) 

where (K 2 ) AB = g CD K AC K BD and TrK 2 = g AB g CD K AC K BD = H 2 - R is expressed in 
terms of the Ricci curvature of the surface by the Gauss-Codazzi relations. Then, the 
extrinsic mean curvature satisfies the parabolic equation 

° = g AB V A V B H + (TrK 2 )H (7.8) 
at 
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generalizing equation (13.131) to two-dimensional branes. 

Clearly, the structure of these equations is much more complicated than those for the 
evolution of planar curves. Although several general results have been obtained in the 
literature so far, the level of our current understanding is by no means complete. Some 
aspects will be discussed here while considering the dimensional reduction of the mean 
curvature flow in R 3 for special classes of surfaces. 

7.2 Dimensional reduction of the curvature flow 

We present two general classes of surfaces, associated to particular ansatz for their em- 
bedding functions, which allow for consistent reduction of the problem to an effective 
curve shortening problem on the plane. They correspond to cylindrical surfaces and sur- 
faces of revolution as brane models. Specific solutions of different topologies will also be 
discussed in the sequel. It should be mentioned, however, that there are other examples 
of surfaces, like the class of ruled surfaces, which do not admit consistent truncation of 
the mean curvature flow away from fixed points (in that case the helicoid). 

(i) Cylindrical surfaces: The simplest possibility arises for surfaces with embedding 
equations 

X = x(s; t) , Y = y{s; t) , Z = u , (7.9) 

where (x(s),y(s)) is the parametric form of a planar curve C that evolves in time and u 
is the second parameter on the surface. These are cylindrical surfaces of the form Rx C 
and it is entirely obvious that their mean curvature flow in R 3 is equivalent to the mean 
curvature flow of the curve C on the plane perpendicular to the Z-axis. As such, they 
provide a trivial dimensional reduction of the mean curvature flow to lower dimensions. 
Other ansatz may also reduce the problem to lower dimension, as deformations of planar 
curves, but the effective dynamics differs from the ordinary mean curvature flow in R 2 , 
as will be seen shortly. 

Any solution of the mean curvature flow on the plane is elevated to a deforming 
cylindrical brane in R 3 , and vice- versa. Thus, for example, all self-shrinking solutions of 
this type are classified by the cylinder R x S 1 and the self-intersecting surfaces R x T p q 
given in terms of the Abresch-Langer curves T p>q . Also, the self-expanding solution, 
associated to the decay of a wedge on the plane, is elevated to an open surface going out 
of the intersection of two planes in R 3 , meeting on the Z-axis. In the context of quantum 
field theory one has D2-branes in the conformal field theory of three free bosons but the 
third boson essentially acts as spectator in the boundary flow equations. 

(ii) Surfaces of revolution: Next, we consider surfaces of revolution in R 3 described 
in all generality by the embedding equation 

X = y( S ;t) cosO , Y = y(s; t) sm9 , Z = x(s;t) . (7.10) 

They are formed by rigid rotation of a planar curve (x(s),y(s)) around the Z-axis by an 
angle 9 that ranges from to 2tt. The time evolution is encoded into the revolving planar 



60 



curve and the mean curvature flow in R 3 is consistently reduced to planar deformations 
of a certain kind. 

The reduced flow is derived by first computing the tangent vectors to the surface, 

(/ Qy \ / Qy \ Qx \ 

I— Jcostf, (— jsinfl, — 1, e e = {-y sintf, y cos6, 0) , (7.11) 

which give rise to the induced metric with components 

^=(|) 2 +(|) 2 > to = f, 9s e = 0. (7.12) 

The normal vector inward to the brane is 

n = , = (cos#, sin6>, — (f'(x)) (7-13) 

using the notation y = (f(x) for the graph of the underlying planar curve. Then, the 
mean curvature of the surface turns out to be 

H = g AB (V A e B ) ■ h = . 1 ( - ^ - -L-) . (7.14) 

The mean curvature flow in R 3 is dimensionally reduced to the following deformation 
on the plane 

&P_ <p"(x) ]_ r71 ^ 

at i + v ,2 {x) (p(x) ' { ' ' 

which clearly differs from the usual mean curvature flow in R 2 by the extra term l/ip(x). 
This difference is attributed to the extrinsic curvature of the S 1 direction following the 
revolution around the Z-axis. Reparametrizations generated by a vector field £ can also 
be added along the flow, as usual. Note, however, that there is no simple variant of 
equation (13.131) satisfied by the curvature of the underlying planar curve that could be 
used further, as in sections 3 and 4. 

Static solutions are characterized by the equation <p(x)<p''(x) = 1 + (p' 2 {x) and they 
correspond to minimal surfaces of revolution in R 3 . In particular, one obtains 

<f(x) = coshx , (7-16) 

which is the graph of the catenary curve on the plane. The solution is equivalently 
described in terms of Liouville equation 

f"(x) + e fix) = (7.17) 

using the relation 

ip 2 (x) cosh x 
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Then, according to the embedding equations in R 3 , the complete surface is described by 
the algebraic equation 

X 2 + Y 2 = cosh 2 Z (7.19) 

so that the two principal curvatures cancel each other at all points. The catenoid surface 
approximates the one-sheeted hyperboloid X 2 + Y 2 — Z 2 = 1 only for small Z\ it yields 
the well known shape of soap bubbles extended between two parallel circular boundaries. 
Another minimal surface in R 3 is the plane, but it appears in a somewhat singular way 
in the present formalism, as surface of revolution of the planar line x = around the 
Z-axis perpendicular to it. 

Self-similar solutions provide the simplest examples of immersed surfaces that evolve 
by overall time scaling so that the functions (I7.10p have common factorized dependence 
by \^2ct. Their position vector f— (X, Y, Z) satisfies by definition the special relation 

Hn = cf . (7.20) 

For surfaces of revolution it yields the effective planar curve equation for y = ip(x) 



ip"(x) 1 

c(y - V x = - f l /2 - — - , 7.21 

1 + if' [x) (p(x) 



since (x(s;t),y(s;t)) also evolves by overall scaling. Here, we only discuss examples of 
self-shrinkers with c < so that t runs from — oo to and the whole surface shrinks to 
the origin by dilations. We will examine solutions with cylindrical, spherical and toroidal 
topology and refer briefly to some of their consequences. 

In all cases, these surfaces correspond to stationary points of Huisken's functional 
cr2(s)/2 ^det^ dsdO = f e <- 2 (s)+y 2 (s))/2 y ^ ^ dx / ds f + (dy/dsf dsd9 (7.22) 



e 

Af v J AT 

for the normalized mean curvature flow in R 3 . Integration over 9 is performed trivially 
and one is left with an integral over s representing the length of a planar curve with 
appropriate metric. In this context, self-similar solutions are effectively described by 
geodesies in the upper half-plane (x, y), with y > 0, that comes equipped with the metric 

ds 2 = y 2 e< x2+y2 \dx 2 + dy 2 ) . (7.23) 

If the factor y 2 were missing the answer would be the same as for the geodesic inter- 
pretation of scaling solutions on the plane found in section 4.4. The presence of this 
additional factor accounts for the extra term l/tp(x) in equation ( 17. 15ft . differentiating 
the dimensionally reduced equation from the ordinary planar mean curvature flow. 

The first example is provided by self-shrinking cylinders of radius ay/2ct with a = 
l/V - c i which are common to the classes (i) and (ii). In the present context, they 
correspond to solutions of equation (17.211) with (p(x) = a. Next, there is the example of 



self- shrinking spheres of radius aylct with a = y— 2/c. They correspond to solutions 
of equation ( 17.211) with ip(x) = \J o? — x 2 , which represents a semi-circle in (x,y) plane. 
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Comparison between the two solutions shows that spheres shrink faster that cylinders of 
equal initial radius. This is also expected on intuitive grounds since spheres are more 
curved than cylinders of equal radii. Finally, there are self- shrinking doughnuts in R 3 
whose existence was first established in Ref. [HZj ; for a discussion see also Ref. [18J. 
They correspond to simple closed geodesic in the upper half-plane equipped with the 
metric (17.231) . which is symmetric with respect to reflection in the y-axis. The proof 
relies on the so called shooting method and proceeds in several steps that are omitted 
here. Unfortunately, there is no closed formula that describes the corresponding planar 
curve that accounts for such solution. Certainly, it can not be a round circle for this does 
not provide solution to the reduced flow (17.211) . 

The classification of self-similar solutions in R 3 and the formation of singularities 
under the flow are not fully explored in all generality. Apart from the obvious scaling 
solutions R x S 1 , R x T p q , the round S 2 and the self-shrinking doughnut there can be 
many more surfaces of various topologies that may also admit self-intersections. The 
general situation is well understood only for surfaces of positive mean curvature, since 
the sphere is the only compact surface of this kind that evolves by scaling, [48] . There are 
other important differences with the mean curvature flow in R 2 that complicate things 
further. Closed curves embedded in the plane always shrink to a point irrespective of 
their initial shape. Even if the curve is not convex at some initial time, it will become 
convex at later times, jTf) . and then approach the homothetic collapse of the round 
circle towards a point, [46J. This property does not generalize to higher dimensions as 
singularities can arise before the surface has the chance to become convex. 

The existence of self-similar shrinking doughnuts can be employed to provide a qual- 
itative proof of this behavior, [67] . For it suffices to consider a dumbbell in R 2 consisting 
of two large approximately round spheres connected with a thin long cylinder as ini- 
tial configuration. By considering a small self-similar shrinking doughnut that encircles 
the neck of the surface, one easily sees that the doughnut, and hence the neck of the 
dumbbell, will become singular well before the two spherical regions have a chance to 
collapse. 

Similar constructions and arguments apply to hypersurfaces of codimension 1 embed- 
ded in all higher dimensional flat spaces. Much less is known about the general features 
of evolving hypersurfaces in flat space when their codimension is bigger than 1 and/or 
when the ambient space is curved. 
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8 Conclusions 



The quantum field theory of two-dimensional sigma models provides a natural framework 
for the realization of both intrinsic and extrinsic curvature flows. These theories have 
all the necessary geometric ingredients to define the flows. Classically, the target space 
fields as well as the embedding equations for branes are fixed once and for all, but, 
in the quantum theory, they are regarded as generalized couplings that depend on the 
energy scale. Thus, the renormalization group equations of the sigma models induce 
flows that can be computed perturbatively. The first order corrections in a' expansion 
are given by the curvature (intrinsic or extrinsic) and the resulting equations combine 
into a coupled system of Ricci and mean curvature flows. There can be additional fields, 
such as anti-symmetric tensor, dilaton and gauge fields, whose beta functions combine 
with the others into a larger system of flows. The deformations of the bulk couplings 
form a closed system, which is independent of the existence of branes and can be studied 
separately. On the other hand, the deformations of the boundary couplings depend on 
the background in which branes are embedded. The resulting picture puts the boundary 
renormalization group equations on firm mathematical base, as for the bulk equations. 
It also suggests generalizations of the combined Ricci and mean curvature flows in the 
presence of fluxes, via the Dirac-Born-Infeld action, which demand further attention. 

Fixed point configurations are reached when the quantum field theory is conformal. It 
is possible, however, to have non-conformal boundary conditions for branes that deform 
in a conformally invariant background. Then, in this context, ordinary D-branes are char- 
acterized by conformally invariant boundary conditions in a conformal field theory. We 
have examined several interesting examples of either kind in two- and three-dimensional 
ambient spaces using appropriate mini-superspace reductions of the more general prob- 
lem. Even in the simplest case of the conformal quantum field theory of two free bosons, 
represented by the plane, the possibilities for brane evolution are enormous and there is 
no systematic way to solve the associated curve shortening problem in all generality. It 
will be useful to develop new algebraic techniques, as for the Ricci flow in two dimensions, 
which will enable to cast the mean curvature flow into zero curvature form. In the same 
spirit, it will be interesting to investigate all integrable perturbations of a given fixed 
point solution, such as the hair-pin, and associate renormalization group trajectories to 
new infra-red fixed points, in analogy with the integrable perturbations of bulk conformal 
field theories. 

The construction of entropy functionals and their physical interpretation in terms 
of the underlying quantum field theory of Dirichlet sigma models are other directions 
of future research. In principle, one should be able to generalize Huisken's functional, 
mentioned in section 3, to branes deforming in curved ambient spaces with or without 
fluxes. Even the simplest cases corresponding to the target space of exact conformal field 
theories, such as the two-dimensional Euclidean black hole or S 3 stabilized by fluxes, have 
not been considered to this day. It is also natural to expect that the critical points of such 
generalized entropy functionals will help to characterize the singularities of collapsing 
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branes in curved spaces, in analogy with the self-shrinking solutions in flat space. When 
the branes deform in running backgrounds the problem becomes even more interesting 
for there can be branes that become singular before or simultaneously with the metric. 
None of these possibilities have been analyzed before in the mathematics literature and 
the corresponding entropy functionals are yet to be found. It should also be noted in 
this context that the ^-function of boundary flows, [HE], [69], is still awaiting its proper 
mathematical place in the framework of mean curvature flows, as for the c-function of 
bulk flows expressed by Perelman's entropy of Ricci flows, [13]. 

The boundary state formalism of Dirichlet sigma models should be developed further 
in order to provide exact characterization of the fixed points as well as the running 
solutions of the flow from the world-sheet view-point. In this context, it will be interesting 
to consider the effect of instantons on the mean curvature flow, as for the Ricci flow, and 
investigate the emergence of non-trivial infra-red fixed points. A simple example of this 
kind is the 0(3) sigma model with 9 = n topological term in which there can be embedded 
closed curves that normally deform to a point. The #-term yields the Gaussian model 
of a free boson as infra-red limit of the bulk theory, which is compactified on a circle of 
self-dual radius, and the branes ought to flow to D-branes on this circle (see, for instance, 
[TO] for their complete classification). Another class of models is provided by the planar 
Abresch-Langer curves, which, in the presence of the appropriate #-term, may give rise 
to some kind of minimal (p, q) exact boundary states in the quantum field theory of two 
free bosons. Other interesting applications arise in the context of tachyon condensation 
in string theory and in the Kondo effect (for a recent discussion see, for instance, |71j). 
where boundary renormalization group equations play pivotal role. All these questions 
are currently under investigation and further results will be reported elsewhere. 



Acknowledgements 

This work was supported in part by the European Research and Training Network 
"Constituents, Fundamental Forces and Symmetries of the Universe" under contract 
number MRTN-CT-2004-005104, the INTAS programme "Strings, Branes and Higher 
Spin Fields" under contract number 03-51-6346, and the EUAN programme of the Gen- 
eral Secretariat for Research and Technology of Greece under contract number B.545. 
C.S. also acknowledges partial support from the programme "Particle Physics Phe- 
nomenology, NCSR-D". I.B. is particularly thankful to the participants of the Workshop 
on Geometric and Renormalization Group Flows, held in Golm, Germany, for many useful 
discussions and to Gerhard Huisken for his interest in this work and encouragement. 



65 



APPENDICES 



A Embedding equations in Riemannian geometry 

In this appendix we review the main parts from the theory of embedding hypersurfaces in 
Riemannian geometry. The presentation is kept quite general so that it can accommodate 
branes of Dirichlet models with arbitrary codimension defined in general ambient target 
spaces. A more complete account can be found in the textbooks; see, for example, the 
classic reference |72j . 

Consider a Riemannian manifold M of dimension m with local coordinate system 
and metric G^X) so that its line element is 

ds 2 M =G iiV {X)dX»dX v ; (jt,i/ = l,2,---,m . (A.1) 

Also consider a submanifold M of M. with dimension n < m and local coordinates y A 
that describes an embedded hypersurface with defining relations X M = f^(y A ). The line 
element in M is given by the corresponding metric gAB{y), 

ds 2 N = g AB (y)dy A dy B ; A, B = 1, 2, • • • , n , (A.2) 

which, of course, is obtained by restricting the line element of the ambient space to N '. 
Thus, gAB^y) is the induced metric on Af equal to 

9ab(v) = G,J^ B . (A.3) 

The tangent vectors to the hypersurface are given in terms of the derivatives of the 
embedding functions, df^/dy A = and they are n of them labeled by the index A. 
Since f 1 are scalars with respect to covariant differentiation on A/", we have equivalently 
DaJ 11 = f^A ■ The (unit) normal vectors to the hypersurface will be denoted by n£, thus 
being labeled with the index a = n + 1, n + 2, • • • , m, and they are chosen to satisfy the 
orthonormalization conditions 

G^X = Sot ■ (A.4) 

By definition they are orthogonal to the tangent vectors to the hypersurface, i.e., 

GWX = . ( A - 5 ) 
and all together they satisfy the following completeness relation in A4, 

9 AB f" A f: B + KKS° T = . (A.6) 

Apart from the induced metric qab there is also the second fundamental (quadratic) 
form on A/", which is a collection of symmetric tensors defined as 

K a AB = GV< (D A D B r + T;j p A f x B ) (A.7) 
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and labeled by the number of transverse directions to the hypersurface. The eigen- values 
of the matrix representing the second fundamental form provide the principal curvatures 
of the hypersurface at each point. The (extrinsic) mean curvature of M in M. is defined 
by taking the trace of the second fundamental form, 

H° = g AB K\ B , (A.8) 

whereas the mean curvature vector associated to each point of the hypersurface is defined 
to be H a n%. When the codimension of the hypersurface is bigger than one there is also 
the so called third fundamental form on A/", which is defined as 

TT = G, U K [h^A + ^Ja) (A.9) 

and it is anti-symmetric under the interchange of its two labels a and r. 

According to these definitions, the equations for the derivatives of the tangent and 
unit normal vector fields on M take the form 

D A D B f» = K\ B K - V^ A f% , (A. 10) 

and 

<a = ^AB9 BC f% - r%fffi - TTK , (A.ll) 

thus extending the Serret-Frenet relations for embedded curves in R 3 to general situa- 
tions. They are called embedding equations since the hypersurface is completely specified 
by the set of these vectors. 

Finally, there are compatibility conditions for the existence of solutions to the em- 
bedding equations, when a given system of tangent and unit vectors is prescribed, which 



take the following form: 

RABCD = + K C[A K B]D > ( A - 12 ) 

D [C K B]A = R» vpX f» A f%fiK + T{SK B]A , (A. 13) 

D [B TXJ + T$T% + g CD KZ [B K A]D + R, vpX f^ B Kfi x T = . (A.14) 



Summation is implicitly assumed over all repeated indices. The first two conditions are 
known as Gauss-Codazzi equations, whereas the last one is known as Ricci equation. They 
all relate the various fundamental forms with the Riemann curvature tensor on M and M.. 
In general, there are more unknown functions than relations in the embedding equations, 
but their number is reduced by performing local transformations in the normal space 
to the hypersurface, which rotate the components of the second and third fundamental 
forms. 
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B Deforming curves and integrability 



The mean curvature flow on the plane is special case of more general dynamics of curves 
that deform as 

dr 

= Un + Wt. (B.l) 

at 

The unit normal and tangent vectors provide the orthonormal base to decompose vectors 
at each point of the curve and the coefficient functions U and W are taken to be local 
functionals of S = — r ■ n or functionals of the extrinsic curvature H and their derivatives. 
It is always convenient to think of the extrinsic curvature as function of the slope /3, in 
which case H((3) is related to S((3) by equation ( 13. lip . Then, the mean curvature flow, in 
its simplest form, df/dt = Hh, corresponds to the choice U = H and W = 0, whereas the 
normalized mean curvature flow df/dt = Hh + f corresponds to U = H(/3) — S(/3) and 
W = S'(P), since f = (f-t)t+ (f-n)n = S'(j3)t — S({3)n. Under the general circumstances 
(IB.ip . the evolution for the extrinsic curvature becomes, [73J, 

^ = H 2 ([/"(/?) +U) + HH'(P) (W"((3) + W) , (B.2) 

substituting for equation (13.131) or (13. 14|) when U(f3) and W(f3) are arbitrary. 

Neither variant of the mean curvature flow preserves the length of the curve. If we are 
prepared to study generalized evolution equations ( IB. If) that keep invariant the length 
of the curves, not only globally but also locally, then we are led to consider coefficient 
functions that satisfy the special relation 

W'{p) = U((3) . (B.3) 

In this case, the derivatives with respect to the arc-length I and time t commute. Fur- 
thermore, if U is of the general form HV'(f3) = dV(l)/dl for some local functional V, the 
evolutions so defined will also preserve the total area surrounded by such closed curves. 
The physical picture is to consider inextensible strings, open or closed, that deform on 
the plane by those general rules. It is quite remarkable that there is an infinite hierarchy 
of flows, other than the mean curvature flow, which satisfy the constraints mentioned 
above and yield integrable equations for H(/3), which in turn can determine the curve 
by equation (I3.12p . This puts our investigation in a wider framework and points out 
that, unlike other cases, the mean curvature flow does not seem to reduce to some known 
integrable system. 

Elaborating; more on this point, note that the choice W(J3) = H 2 (/3)/2 and U(J3) = 
HH'([3), which is consistent with the relation (1B.3I) . amounts to having a third order 
differential equation for H(/3), as follows from the general evolution (1B.2I) above, that 
is equivalent to the modified Korteweg-de Vries (mKdV) equation. This is best seen in 
terms of the arc-length I of the curve, rather than /3, since the corresponding equation 
for the extrinsic curvature takes the standard mKdV form 



^ = H"\l) + -H*H'(l) , (B.4) 
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where prime denotes here the derivative with respect to I. In this context, W is naturally 
identified with the first conserved quantity of the mKdV equation. More generally, one 
can show that choosing W as any one of the higher conserved quatities of the mKdV 
equation (1B.4j) . and U according to the relation (IB.3I) . amounts to reducing the general 
evolution (1B.2|) for the curvature into integrable equations that coincide with the other 
members of mKdV hierarchy. 

Finally, we note for completeness that several other type of integrable systems have 
been obtained in the literature by considering evolutions of curves in two or higher 
dimensional spaces, which may also be curved, via general evolutions of the form ( IB. II) . 
None of these, however, can accommodate the mean curvature flow in two or higher 
dimensional ambient spaces. 



C Resistive diffusion of magnetic fields 

In this appendix we review the emergence of the mean curvature flow for planar curves 
via dimensional reduction of the magneto-hydrodynamic equations for time dependent 
force-free magnetic fields in R 3 , following Ref. |54j . 

Consider the resistive diffusion of a magnetic field B(x, t) in a medium (plasma), which 
is free to move with velocity v(x, t) to accommodate the changing magnetic configuration 
in real time. The basic magneto-hydrodynamic equation controlling the process is 

-+i)Vx(Vx5) = Vx(t?x5), (C.l) 
at 

where rj is the constant resistivity of the medium. It provides a good approximation to 
the real world in the limit of vanishing gas pressure when the magnetic field obeys the 
force-free condition (V x B) x B = 0. Thus, one is led to consider magnetic fields of the 
form 

VxB = aB, (C.2) 
where a(x,t) is a scalar function satisfying the special relation 

(B-V)a = (C.3) 

so that the source free Maxwell equation V • B = is obeyed. The case of constant a in a 
static medium is trivial and corresponds to an exponentially decaying field satisfying the 
linear equation dB/dt + r)a 2 B = 0. In the following we will be concerned with magnetic 
fields with non-constant a. 

The relations fIC.ll) . (10. 2[) and (10.31) form, in general, a non-linear coupled system of 
seven equations for seven unknown (B,v,a) that fully determine the resistive diffusion 
of a force-free magnetic field in passive medium. Here we consider their reduction to one 
spatial direction, say z, by assuming that the magnetic field takes the special form 

B(z,t) = {B cosl3(z,t), B smf3{z,t), 0), (C.4) 
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where B is constant that can be set equal to 1, whereas the velocity in the medium is 
taken to be 

i?=(0, 0, v z (z,t)). (C.5) 

Then, the magneto-hydrodynamic equations reduce to the following simpler system for 
the unknown functions (3 and v z , 

whereas the function a, which is also considered as function of z and t, is determined by 

a(z,t) = -^. (C.8) 

Thus, one is only left to determine (3 and v z . 

Considering v z (z,t) as function of f3(z,t) and t, and defining the quantity 

ff(At) = ^p, (C9) 



it follows that the system of equations (1C.6j) and (1C.7j) take the equivalent form 



V ^ + H((3,t) = 0, (CIO) 

According to the ansatz (10.41) . the magnetic field B(x,t) is uniform on the z = constant 
planes and changes direction as one moves across these planes. Close inspection with 
the mean curvature flow of planar curves shows that —z/rj represents the arc-length of a 
virtual curve, as measured from the origin of coordinates, z = 0, j3 is the slope and H(/3) 
the extrinsic curvature at each point of the curve. The magnetic field on the z = constant 
planes, B Q (cos/3, sin/3), provides the tangent vector to such curves at different points and 
the function a(z) equals to H/r]. Finally, the velocity of the medium, v z (z), is nothing 
else but the (Euler-Bernoulli) elastic energy of the curve segment of length z/rj, assuming 
that the velocity vanishes at z = 0. Then, in this context, the evolution equation (1C.11I) 
coincides with equation ( 13.131) for the extrinsic curvature of the deforming virtual curve 
with respect to the rescaled time t/rj. Any solution of the mean curvature flow gives rise 
to a process for the resistive diffusion of force-free magnetic fields in R 3 . 

A typical initial value problem for these equations is defined by the following condi- 
tions 

/3{z, 0) = -p{-z, 0) , ^j^>0 (C.12) 
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for t = 0, as — oo < z < oo, together with the boundary conditions 

P(±oa,t) = ±P (C.13) 

for all t > 0, having finite (3q independent of t; the initial and boundary conditions for 
v z (z, t) follow from equation (1C6|) . The direction of B is held fixed at z = ±00 for all time 
by the above boundary conditions, experiencing an overall rotation by 2/?o from one end 
to the other. An interesting question that arises in this context, and has applications in 
modeling the eruption of solar flares, is to determine the conditions under which H(j3, t) 
becomes infinite at some time. According to the analysis of Ref . [54J , a modest criterion 
is provided by comparing 2/5 to ir. If 2f3 Q < ir, B(z,t) will evolve toward a uniform 
field irrespective of the initial value /3(z,0). If 2/3q > tt, B(z,t) will develop infinite field 
gradient after some time, irrespective of the initial value /3(z,0), thus leading to infinite 
H by equation ( 1C.10I) . The marginal case 2/3 = vr includes the steady state solution 

B(z) = {±B sech(A2), B tanh(Az), 0) , (C.14) 

v z (z) = -r] A tanh(Az) , a(z) = sech(Az) (C.15) 
that corresponds to the translating soliton of the mean curvature flow on the plane. 

It is quite interesting, in many respects, that apart from the translating soliton other 
simple solutions of these equations were constructed in Ref. [51], including the paper-clip 
model (in modern language). Scaling solutions (self-shrinkers and expanders) were also 
investigated there to some extend. There is appropriate mention to them in the main 
text of our paper. 
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